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1 Introduction 



At the beginning of 2011 the status of the inclusive cross-section for Higgs-boson production in gluon 
fusion was summarized [T]. Corrections arising from higher-order QCD, eleetroweak effects, as weh as 
contributions beyond the commonly-used effective theory approximation were analyzed. Uncertainties arising 
from missing terms in the perturbative expansion, as well as imprecise knowledge of parton distribution 
functions, were estimated to range from approximately 15—20%, with the precise value depending on the 
Higgs boson mass. For an updated study we we refer to Ref. 2 . 

In this work we consider the problem of a consistent definition of the Higgs-boson-lineshape at LHC and 
address the question of optimal presentation of LHC data. 

Some preliminary observations are needed: the Higgs system of the Standard Model has almost a non- 
perturbative behavior, even at a low scale. Let us consider the traditional on-shell approach and the values 
for the on-shell decay width of the Higgs boson [T]: for = 140 GeV the total width is 8.12 x 10"^ GeV; 
the process H — s- bb has a partial with of 2.55 x 10~^ GeV, other two-body decays are almost negligible while 
H — >■ 4f has 4.64 x 10^^ GeV well below the WW -threshold. Therefore the four-body decay, even below 
threshold, is more important than the two-body ones. What are the corresponding implications? Since decay 
widths are related to the imaginary parts of loop diagrams the above statement is equivalent to say that, in 
terms of imaginary parts of the H self-energy, three-loop diagrams are as important as one-loop diagrams, 
or more. 

We realize that most of the people just want to use some well-defined recipe without having to dig any 
deeper; however, there is no alternative to a complete description of LHC processes which has to include the 
complete matrix elements for all relevant processes; splitting the whole 5' -matrix element into components 
is just conventional wisdom. However, the precise tone and degree of formality must be dictated by gauge 
invariance. 

The framework discussed in this paper is general enough and can be used for all processes and for all 
kinematical regions; however, the main focus of this work will be on the description of the Higgs-boson- 
lineshape in the heavy Higgs region, typically above 600 GeV. The general argument is well documented 
in the literature and, for a complete description of all technical details, we refer to Ref. |3I4| . Part of the 
results discussed in this paper and concerning a correct implementation of the Higgs-boson-lineshape have 
been summarized in a talk at the BNL meeting of the HXSWG0. It is important to stress the exact content 
of this work; we deal with a gauge invariant definition of the signal, which assumes that the non-resonant 
background and the signal/background interference have been subtracted from the data. 

One might wonder why considering a Standard Model (SM) Higgs boson in such a high-mass range. 
There are classic constraints on the Higgs boson mass coming from unitarity, triviality and vacuum stability, 
precision eleetroweak data and absence of fine-tuning However, the search for a SM Higgs boson over a 
mass range from 80 GeV to 1 TeV is clearly indicated as a priority in many experimentalpapers, e.g.ATLAS: 
letter of intent for a general-purpose pp experiment at the large hadron collider at CERI^. For recent results, 
see Ref. |6l7l8l9ll0llirT2] . 

The situation is different if we consider extensions of the SM: in the two Higgs doublet (THD) model, 
even if the SM-like Higgs boson is found to be light (< 140 GeV), there is a possible range of mass splitting in 
the heavy Higgs boson. In general, for a given Higgs boson mass, the magnitude of the mass splittings among 
different heavy scalar bosons can be determined to satisfy the eleetroweak precision data, see Ref. 131. 

In addition, the work of Ref. [2] studies scenarios where a heavy Higgs boson can be made consistent 
with both the indirect constraints and the direct null searches by adding only one new particle beyond the 
Standard Model. Heavy Higgs effects in a theory with a fourth Standard-Model-like fermion generation have 
been examined in Refs. [15|16j . 

The work of Ref. [17] critically reviews models that allow a heavy Higgs boson consistent with the 
precision eleetroweak constraints. All have unusual features, and all can be distinguished from the Minimal 
SM either by improved precision measurements or by other signatures accessible to next-generation colliders. 

-"^http: / /www. bnl.gov/hcs/ 

■^http: / / atlas, wob.ccrn.ch / Atlas / internal / tdr .html 
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Coming back to the framework that we are introducing, there is another important issue: when working 
in the on-sheh scheme one finds that the two-loop corrections to the on-sheU Higgs width exceed the one-loop 
corrections if the on-shell Higgs mass is larger than 900 GeV, as discussed in Ref. iTBl. This fact simply 
tells you that perturbation theory diverges badly, starting from approximately 1 TeV. In this work we 
will also illustrate the corresponding impact on the Higgs boson lineshape (previous work can be found in 
Refs. [I9I20] V 

Recently the problem of going beyond the zero-width approximation has received new boost from the 
work of Refs. }21|22j : the program iHixs allow the study of the Higgs-boson- lineshape for a finite width 
of the Higgs boson and computes the cross-section sampling over a Breit-Wigner distribution. There is, 
however, a point that has been ignored in all calculations performed so far: the Higgs boson is an unstable 
particle and should be removed from the in/out bases in the Hilbert space, without destroying the unitarity 
of the theory. Therefore, concepts as the production of an unstable particle or its partial decay widths do not 
have a precise meaning and should be replaced by a conventionalized definition which respects first principles 
of Quantum Field Theory (QFT). 

This paper is organized as follows. In Section [2] we introduce and discuss complex poles for unstable 
particles. In Section |3] we analyze production and decay of a Higgs boson at LHC. A discussion on gauge 
invariance is presented in Section |4l In Section [6] we present a short discussion on the QCD scale error. 
In Section [5] we present numerical results while in Section [7] we discuss the residual theoretical uncertainty. 
Finally, technical details are discussed in Appendices, in particular in Appendix [B] we discuss how to apply 
the equivalence theorem for virtual vector-bosons and in Appendix I C . II we discuss analytic continuation in 
a theory with unstable particles. 



2 Propagation 

To start our discussion we consider the process ij — )■ H(— > F) -|-X where i,j € partons and F is a generic 
final state (e.g. F = yy, 4f, etc.). For the sake of simplicity we neglect, for a moment, folding the partonic 
process with parton distribution functions (PDFs). Since the Higgs boson is a scalar resonance we can split 
the whole process into three parts, production, propagation and decay. In QFT all amplitudes are made out 
of propagators and vertices and the (Dyson-resummed) propagator for the Higgs boson reads as follows: 

(1) 

where Mi is a renormalized mass and 5hh is the rcnormalized Higgs self-energy (to all orders but with 
one-particle-irreducible diagrams) . The first argument of the self-energy in Eq. ([IJ is the external momentum 
squared, the remaining ones are (renormalized) masses in the loops. We define complex poles for unstable 
particles as the (complex) solutions of the following system: 

SH - Ml + 5hh (sh, Af2, Ml,M^,Mi^ = 0, 
sw - M^ + Sww [sw,Ml Ml,M^,Mi^ = 0, (2) 
etc. To lowest order accuracy the Higgs propagator can be rewritten as 

Ah'=s-sh. (3) 
The complex pole describing an unstable particle is conventionally parametrized as 

St= -iHiji, (4) 

where /ii is an input parameter (similar to the on-shell mass) while 7^ can be computed (as the on-shell 
total width), say within the Standard Model. There are other, equivalent, parametrizations [18^, e.g. y^sh — 
fin — «/2 7h. Note that the the pole of A fully embodies the propagation properties of a particle. We know 



Ah(s) 



Shh f s. 



M,^M^,M^,Mi 
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that even in ordinary Quantum Mechanics the resonances are described as the complex energy poles in the 
scattering amplitude. The general formalism for describing unstable particles in QFT was developed long ago, 
see Refs. |23I24I25I26I27| : for an implementation in gauge theories we refer to the work of Refs. [28129130131] . 
for complex poles in general to Refs. [32133134] and Refs. [35136] . For alternative approaches we refer to the 
work of Refs. [37l38] . 

We can summarize by saying that unstable particles are described by irreducible, nonunitary representa- 
tions of the Poincare group, corresponding to the complex eigenvalues of the four-momentum satisfying 
the condition = —fi^ + i /i7. For a discussion on violation of the spectral condition we refer to Ref. [39 . 

Consider the complex mass scheme (CMS) introduced in Ref. [40141] (see also Ref. [15]) and extended at 
two- loop level in Ref. [43]: here, at lowest level of accuracy, we use 



M2 



Shh (sh, St, sh, sw, sz) 



0, 



(5) 



where now S'hh is computed at one loop level and sw stc. are the experimental complex poles, derived from 
the corresponding (experimental) on-shcU masses and widths. For the W and Z bosons the input parameter 
set (IPS) is defined in terms of pseudo-observables; at first, on-shell pseudo-quantities are derived by fitting 
the experimental lineshapes with 



(s) = 



(s 



= W,Z, 



(6) 



where N is an irrelevant (for our purposes) normalization constant. Secondly, we define pseudo-observables 



Mp — Afos cos-^, 



Fp = r 



OS 



sinV', 



ip — arctan 



F 



OS 



M, 



OS 



(7) 



which are inserted in the IPS. Renormalization with complex poles should not be confused with a simple 
recipe for the replacement of running widths with constant widths; there are many more ingredients in the 
scheme. It is worth noting that perturbation theory based on sh instead of the on-shell mass has much 
better convergence properties; indeed, as shown in Ref. [18] . the two- loop corrections to the imaginary part 
of Sh become as large as the one-loop ones only when /xh = 1.74 TeV. This suggests that the complex pole 
scheme is preferable also from the point of view of describing the heavy Higgs-boson production at LHC. 

There is a substantial difference between W, Z complex poles and the Higgs complex pole. In the first 
case W, Z bosons decay predominantly into two (massless) fermions while for the Higgs boson below the 
WW -threshold the decay into four fermions is even larger than the decay into a bb pair, while the other 
fermion pairs are heavily suppressed by vanishing Yukawa couplings. Therefore we cannot use for the Higgs 
boson the well-known result. 



Im S'vv(s) 



pOS 
^ V 



which is valid for W, Z. As a consequence of this fact we have 



(8) 



OS 



vos 



HO, 7v = F 



OS 
V 



2 U/OS/ 



HO, 



(9) 



where the perturbative expansion is well under control since Fy^/My^ ^ 1. For the Higgs boson we have a 
different expansion. 



ao = - [Im 5'hh , os] ^ ^ ['^ ^™ "^hh , os " ^^h Im S'hh , os 



"S'hh , OS + S'hh , os ^hh , os I™ ^hh , os ^hh , os 
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•^HH^'os = ...QM^ [Mi,Mt, MA, Mi,, Mi j , (11) 



&n = ^ I™ ^HH . OS - Mh Im 5^^^ , OS , = " ^ ^ Im ^Sh , os (10) 

Here we have Mi = M^^, Ti = Tf^ (on-shell quantities) and 

d 

where only the first few terms in the perturbative expansion are exphcitly given. In Eq. ljlOp we used the 
following definition of the on-shell width 

rg^A^H^-Im^HH,os+HO (12) 

Only the complex pole is gauge-parameter independent to all orders of perturbation theory while on-shell 
quantities are ill-defined beyond lowest order. Indeed, in the gauge, at lowest order, one has the following 
expression for the bosonic part of the Higgs self-energy: 



Im5HH.bos(s) = -^sM^-l l-4ew— .s-4ewM^ +- (W^Z) , (13) 



where (V ~ W, Z) are gauge parameters. Note that Ea.([TUl) (the expansion) involves derivatives. 

Starting from Eq. (jl3p we understand the problem with the on-shell definition of the mass: once again, 
let Mh be the renormalized Higgs mass and sh = /^h ~ * Mh7h the parametrization of the complex pole; it 
follows 

Ml |os= Ml + Re 5hh(m2), ^\ ^ Ml + Re 5hh(sh), 

Mh Ios -m = - ^7H Im S'hh(mh). (14) 

showing that at order the on-shell mass is ill-defined since sh is ^-independent. 

A technical remark: suppose that we use /in as a free input parameter and derive 7h in the SM from the 
equation 

mm = Im Shh (sH,St,SH,sw,sz) . (15) 

Since the bare W, Z masses are replaced with complex poles also in couplings (to preserve gauge invariance) 
it follows that 7h (solution of Ea. (|15p ) is renormalization scale dependent; while counterterms can be in- 
troduced to make the self-energy ultraviolet finite, the /iR dependence drops only in subtracted quantities, 
i.e. after finite renormalization, something that would require knowledge of the experimental Higgs complex 
pole [44l43j . Following our intuition we fix the scale at fi^ = /in • 

After evaluating the coefficients of Eg. dTUl) (e.g. in the ^ = 1 gauge) it is easily seen that the expanded 
solution of Eq.(l5]) is not a good approximation to exact one, especially for high values of nn. For instance, for 
/iH — 500 GeV we have an exact solution 7h — 58.70 GeV, an expanded one of 62.87 GeV with an on-shcU 
with of 68.0 GeV . Therefore, we will use an exact, numerical, solution for 7h; details on the structure of the 
Higgs self-energy to all orders of perturbation theory can be found in Ref. |44I43| . 

It is important to remark that an evaluation of 7h at the same level of accuracy to which F^^ is known 
(see Ref. |^1 ) would require, at least, a three- loop calculation (the first instance where we have a four-fermion 
cut of the H self-energy). 

Complex poles for unstable W, Z, H and t also tell us that it is very difficult for an heavy Higgs boson to 
come out right within the SM; these complex poles are solutions of a (coupled) system of equations but for 
W, Z and (partially) t we can compare with the corresponding experimental quantities. Results are shown in 
Table[l]and clearly indicate mismatches between predicted and experimental W, t complex poles. Indeed, as 
soon as /ifj increases, it becomes more and more difficult to find complex W, t and Z poles with an imaginary 
part compatible with measurements. 
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Table 1: The Higgs boson complex pole at fixed values of the W, t complex poles compared with the complete solution 
for sh , sw and St 



/iH [GeV] 7w [GeV] fixed 7t [GeV] fixed 7h [GeV] derived 

200 2.088 1.481 1.355 

250 3.865 

300 8.137 

350 14.886 

400 26.598 

/iH [GcV] 7w [GeV] derived 7t [GeV] derived 7h [GeV] derived 

200 2.130 1.085 1.356 

250 2.119 0.962 3.823 

300 2.193 0.836 8.139 

350 2.607 0.711 14.653 

400 3.922 0.566 25.498 



This simple fact raises the following question: whai is the physical meaning of an heavy Higgs boson 
search? We have the usual and well-known considerations [S]: a Higgs boson above 600 GeV requires new 
physics at 1 TeV, argument based on partial- wave unitarity |45I46) (which should not be taken quantitatively 
or too literally); violation of unitarity bound possibly implies the presence of J = 0, 1 resonances but there is 
no way to predict their masses, simply scaling the tz—h system gives resonances in the 1 TeV range. Generally 
speaking, it would be a good idea to address this search as search for J = 0, 1 heavy resonances decaying 
into VV — >■ 4f. In a model independent approach both fin and 7h should be kept free in order to perform a 
2 dim scan of the Higgs-boson lineshape. For the high-mass region this remains our recommended strategy. 
Once the fits are performed it will be left to theorists to struggle with a model-dependent interpretation of 
the results. 

To summarize, we have addressed the following question: what is the common sense definition of mass 
and width of the Higgs boson? We have several options. 



SH=A*H-«m7H, SH = ( A*H - -^Th ) > SH = _2 '_2 " ■ (16) 

We may ask which one is correct, approximate or closer to the experimental peak. Here we have to distinguish: 
for 7h ^ /^H Afn is a good approximation to the on-shell mass and it is closer to the experimental peak; 
for instance, for the Z boson Mz is equivalent to the mass measured at Lep. However, in the high-mass 
scenario, where 7h ^ I-I-h, the situation changes an On ^ for any observable (or pseudo-observable). 

Therefore, the message is: do not use the on-shell width to estimate Afn in the high-mass region. 

Missing a calculation of the three-loop Higgs self-energy we can analyze the low-mass region by using an 
approximation to the exact Higgs complex-pole which is based on the expansion of 7h around the best-known 
on-shell calculation, Tp from Prophecy4f [47|. Therefore, below = 200 GeV we will use 

l^^l + lxw (di-fil,d2) ^-X^dl (17) 

Xw=4V2 d„ = — 5hh. (18) 

In the expansion all masses but the Higgs mass are kept real and the higher-order effects are simulated by 
the expansion parameter, Tp/ fin- 
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3 Production and decay 



Before describing production and decay of an Higgs boson we underline the general structure of any 
process containing a Higgs boson intermediate state. The corresponding amplitude is schematically given by 



S - SH 



(19) 



where N{s) denotes the part of the amplitude which is non-Higgs- resonant. Signal (S) and background (B) 
are defined as follows: 



Ms) = Sis) + Bis), 5(.) = /M, Bis)^l^^t^I^ + Nis) 



S - SH 



(20) 



As a first step we will show how to write /(s) in a way such that pseudo-observables make their appearance. 
Consider the process ij — H — F where i,j G partons and F is a generic final state; the complete cross- 
section will be written as follows: 



1 

2^ 



A. 



1 



S - SH 



(21) 



where ^ is over spin and colors (averaging on the initial state). Note that the background (e.g. gg — 4f, 
see Ref. }48l49l50j ) has not been included and, strictly speaking and for reasons of gauge invariance, one 
should consider only the residue of the Higgs- resonant amplitude at the complex pole, as described in Eg. ipOl) . 
For the moment we will argue that the dominant corrections are the QCD ones where we have no problem 
of gauge parameter dependence. If we decide to keep the Higgs boson off-shell also in the resonant part of 
the amplitude (interference signal/background remains unaddressed) then we can write 



sAijis). 



(22) 



For instance, we have 



Gp s 



TT^ 288 



(23) 



where Tg = Arn^/s, f{Tq) is defined in Eq.(3) of Ref. [21] and where (5qcd gives the QCD corrections to 
gg — >■ H up to next-to-next-to-leading-order (NNLO) + next-to-leading logarithms (NLL) resummation. 
Furthermore, we define 

(24) 

S.C 



rH^F(s) = / d*H^F EI^H^F 

which gives the partial decay width of a Higgs boson of virtuality s into a final state F. 

A,j(s) 



(25) 



which gives the production cross-section of a Higgs boson of virtuality s. We can write the final result in 
terms of pseudo-observables 



O'ij-i.H^F(s) = —au^u 

IT 



It is also convenient to rewrite the result as 



0'i7"^H-i.F(s) = —(^ij^H 
TT 



S - SH 



T^tot 

^ H 



S - SH 



— ^BR (H 



(26) 



(27) 
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where we have introduced a sum over all final states, 



= E ^H^F. (28) 

Note that we have written the phase-space integral for i{pi) + j{p2) -> F as 

J d$,y^F= I d^k6%k-p,-p2) J l[d*pfS+{p})5\k~J2Pf)' (29) 

where we assume that all initial and final states (e.g. yy, 4f, etc.) are massless. 

It is worth noting that the introduction of complex poles does not imply complex kinematics. According 
to Eg. lpO)) only the residue of the propagator at the complex pole becomes complex, not any element of the 
phase-space integral. Details of the procedure are discussed in Appendix [U] 

On a more formal bases one should say that unstable states lie in a natural extension of the usual 
Hilbert space that corresponds to the second sheet of the S* -matrix; these states have zero norm and, 
therefore, escape the usual prohibition of having an hermitian Hamiltonian with complex energy 1521 . On a 
more pragmatic level we use the guiding principle that Green's functions involving unstable particles should 
smoothly approach the value for stable ones (the usual Feynman — i prescription) when the couplings of 
the theory tend to zero. 

If we insist that |H > is an asymptotic state in the Hilbert space then the observable to consider will 
be < F out I H in >, otherwise one should realize that for stable particles the proof of the LSZ reduction 
formulas | 53| depends on the existence of asymptotic states 

Ipin lim [ d^xH{x)iSte'P''\0 >, (30) 

(in the weak operator sense). For unstable particles the energy is complex so that this limit either diverges 
or vanishes. Although a modification of the LSZ reduction formulas has been proposed long ago for unstable 
particles, see Ref. [52] . we prefer an alternative approach where one considers extracting information on the 
Higgs boson directly from 

< F out I H >< H 1 1 in > + E < F out I n >< n 1 1 in >, (31) 

for some initial state I and some final state F and where {n} ® H is a complete set of states (not as in 
the in/out bases). As we have seen, the price to be paid is the necessity of moving into the complex plane. 
Why do we need pseudo-observables? Ideally experimenters (should) extract so-called realistic observables 
from raw data, e.g. a (pp — > yy -I- X) and (should) present results in a form that can be useful for comparing 
them with theoretical predictions, i.e. the results should be transformed into pseudo-observables; during the 
deconvolution procedure one should also account for the interference background - signal; theorists (should) 
compute pseudo-observables using the best available technology and satisfying a list of demands from the 
self-consistency of the underlying theory. 

• Definition We define an off-shell production cross-section (for all channels) as follows: 

1 s2 rtf' 



SH 



When the cross-section ij — >■ H refers to an off-shell Higgs boson the choice of the QCD scales should be 
made according to the virtuality and not to a fixed value. Therefore, always referring to Figure ^ for the 
PDFs and aij^-a+x one should select /i| = /i^ = zs/4 {z s being the invariant mass of the detectable final 
state). Indeed, beyond lowest order (LO) one must not choose the invariant mass of the incoming partons for 
the renormalization and factorization scales, with the factor 1/2 motivated by an improved convergency of 
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fixed order expansion, but an infrared safe quantity fixed from the detectable final state, see Ref. j54] . The 
argument is based on minimization of the universal logarithms (DGLAP) and not the process-dependent 
onefjfl. 

The off-shell Higgs-boson production is currently computed according to the replacement 

^os {^^l)S{C-^il)=^ aoFS (C) BW(C) , (33) 

(e.g. see Ref. |55j) at least at lowest QCD order, where the so-called modified Breit-Wigner distribution is 
defined by 

BW(.)^-- (34) 

where now /in — M^^ and — zs. This ad-hoc Breit-Wigner cannot be derived from QFT and also is 
not normalizable in [0 , +oo]. Note that this Breit-Wigner for a running width comes from the substitution 
of r — r(s) = Fs/M^ in the Breit-Wigner for a fixed width F. This substitution is not justifiable. Its 
practical purpose is to enforce a physical behavior for low virtualities of the Higgs boson but the usage cannot 
be justified or recommended. For instance, if one considers VV scattering and uses this distribution in the 
s -channel Higgs exchange, the behavior for large values of s spoils unitarity cancellation with the contact 
diagram. It is worth noting that the alternative replacement 

^os(mI)5(C-mI)^^ofs(C)BW(C), BW(C) = -- (35) 



" (C-M&) +(A^Hros)" 



has additional problems at the low energy tail of the resonance due to the gg -luminosity, creating an artificial 
increase of the lineshape at low virtualitiesQ. 

Another important issue is that 7h which appears in the imaginary part of the inverse Dyson-resummed 
propagator is not the on-shell width since they differ by higher-order terms and their relations becomes 
non-perturbative when the on-shell width becomes of the same order of the on-shell mass (typically, for 
on-shell masses above 800 GeV). For nonlinear parameterizations of the scalar sector of the standard model 
ande Dyson summation of the Higgs self energy we refer to the work of Ref. ^S5] . 

The complex-mass scheme can be translated into a more familiar language by introducing the Bar-scheme. 
Using Eq.([3]) with the parametrization of Eq.(j4]) we perform the well-known transformation 

= /^H + 7h Mh Fh = A/h 7h. (36) 
A remarkable identity follows (defining the Bar-scheme): 

' =(l + ^i^)Cs-Ml+^I^s)-\ (37) 



showing that the Bar-scheme is equivalent to introducing a running width in the propagator with parameters 
that are not the on-shell ones. Special attention goes to the numerator in Ea. (j37p which is essential in 
providing the right asymptotic behavior when s — >■ cxd, as needed for cancellations with contact terms in 
VV scattering. If we compare the result of Eg. ([571) with Eq.(4) of Ref. [57] and interpret toh and Fh of that 
equation as the corresponding Bar - scheme quantities Mh and Fh of Ea. (p6l) . we see that the so-called 
Seymour - scheme of Refs. j21l22j is exactly giving the Higgs propagator with a complex pole. There is a 
second variant for the Higgs propagator in Ref. [57], i.e. (in the notation of that paper) 



F 



H 



Ah(s) = s-mf,+i s (38) 

TOH 



^We gratefully acknowledge S. Forte and M. Spira for an important discussion on this point. 
^We gratefully acknowledge S. Frixione for an important discussion on this point. 
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Table 2: Higgs boson complex pole; F^^ is the on-shell width, 7h is defined in Eq.Q and the Bar-scheme in Ea. H36p . 



MH[GeV] 


r5J«[GeV] (Ref. [J) 


7H [GeV] 


MH[GeV] 


rH[GeV] 


200 


1.43 


1.35 


200 


1.35 


400 


29.2 


25.60 


400.9 


26.66 


600 


123 


103.93 


608.9 


105.48 


700 


199 


162.97 


718.7 


167.33 


800 


304 


235.57 


834.0 


245.57 


900 


449 


320.55 


955.4 


340.28 


1000 


647 


416.12 


1083.1 


450.71 



which has been used afterwards with the motivation that the prescription accounts for signal-background 
interference effects. If we denote by An the Higgs-resonant amplitude and by the box contribution in 
gg —7- ZZ then the interference is 

Ant - 2 Re (Ah /t) , I = AiiAl, A^ = {s - s^i) A^. (39) 

Note that for large values of the Higgs mass the term in Ea. ([M)) proportional to Im Ah is not negligible. 
The main issue in Ref. is on unitarity cancellations at high energy. Of course, the behavior of both 
amplitudes for s — (X) is known and simple and any correct treatment of perturbation theory (no mixing of 
different orders) will respect the unitarity cancellations; from this point of view, the numerator in Eq.p7|) is 
essential. Since the Higgs boson decays almost completely into longitudinal Zs, ior s oo we have [5^ (for 
a single quark q) 

srai „ s mj: „ s 

showing cancellation in the limit (s Ah — >■ 1). However, the behavior for s — > cx) (unitarity) should not/cannot 
be used to simulate the interference for s < . The only relevant message to be derived here is that unitarity 
requires the interference to be destructive at large values of s, see also Ref. 159; ■ 

A sample of numerical results is shown in Table [5] where we compare the Higgs boson complex pole to 
the corresponding quantities in the Bar-scheme. 

In conclusion, the use of the complex pole is recommended even if the accuracy at which its imaginary 
part can be computed is not of the same quality as the next-to-leading-order (NLO) accuracy of the on- 
shell width. Nevertheless, the use of a solid prediction (from a theoretical point of view) should always be 
preferred to the introduction of ill-defined quantities (lack of gauge invariance) . 

3.1 Schemes 

We are now in a position to give a more detailed description of the strategy behind Eq. (I^Hl) . Consider 
the complete amplitude for a given process, e.g. the one in Figure |9j let {( = zs, . . .) the the full list of 
Mandelstam invariants characterizing the process, then 

A (C, . . .) = Vp,od (C, . . .) Aprop(C) Kicc(C) +N{C,...). (41) 

Here Vprod denotes the amplitude for production, e.g. gg — > H(^) +X, Aprop(C) is the propagation function, 
Vdec(C) is the amplitude for decay, e.g. H(C) — > 4f. If no attempt is made to split A{s) no ambiguity arises 
but, usually, the two components are known at different orders. Ho to define the signal? The following 
schemes are available: 

ONBW 

5 (C, . . .) = Vp,od (a*h, • ■ •) Ap,.op(C) VdU^^l) Aprop(C) = Breit-Wigner. (42) 

in general violates gauge invariance, neglects the Higgs off-shellness and introduces the ad hoc Breit- 
Wigner of Eq.dMl). 
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OFFBW 

5 (C, . . .) = ^prod (C, • • •) Ap,op(C) 14cc(C), Ap,.op(C) = Breit-Wigner, (43) 
in general violates gauge invariance, and introduces an ad hoc Breit-Wigner. 

ONP 

Sic, ■■■) ^ Vprod (mL • • ■) Api.op(C) Vdcc(/i|), Api.op(C) = propagator (44) 
in general violates gauge invariance and neglects the Higgs off-shellness. 
OFFP 

5 (C, . . .) = Vp,od (C, • • •) Aprop(C) VdeciO, Ap„p(C) = propagator, (45) 
in general violates gauge invariance 

CPP 

S {C, ■■■) = Vprod (sH, • • •) Aprop(C) V;jcc(sh), Aprop(C) = propagator. (46) 

This scheme respects all requirements: only the pole, the residue and the reminder of A{C) are gauge 
invariant. Furthermore the CPP-scheme allows to identify POs like the production cross-section and 
any partial decay width by putting in one-to-one correspondence robust theoretical quantities and 
experimental data. 

From the list above it follows that only CPP-scheme should be used, once the signal/background interference 
becomes available. However, the largest part of the available calculations is not yet equipped with Feynman 
integrals on the second Riemann sheet (complex poles lie there); furthermore, the largest corrections in the 
production are from QCD and we could argue that gauge invariance is not an issue over there, so that 
the OFFP-scheme remains, at the moment, the most pragmatic alternative. Additional considerations with 
more detailed descriptions are however postponed until Section |4l 

Implementation of the CPP-scheme requires some care in the presence of four-leg processes (or more). 
The natural choice is to perform analytical continuation from real kinematics to complex invariants; consider, 
for instance, the process g{pi)+g{p2) — >■ g(P3)+H(p4) where the external H leg is continued from real on-shell 
mass to sh- For the three invariants, 

s^-{pi+P2f, t^-ipi-pif, u^-ipi-psf, (47) 
we use the following continuation: 

i--^ (l-^) (1-COS0), « = -f (l-^) (1 + COS0), (48) 

where 9 is the (real) scattering angle (between p2 and ps) and s = 4i?^, where 2 E is the (real) energy 
of center-of-mass system. For the process g(pi) + g(p2) HiPs) + H(p4) where both external H legs are 
continued from real on-shell mass to sh we obtain 



s 

t = - - 
2 



1 (1-4 — I cos( 



s s 



s 

u = 

2 



l-f^i + (l_4f^i)^/%os0l. (49) 

s s 



For a more detailed discussion of the analytic continuation we refer to Appendix IC.ll Here we observe that 
when the problem is with extracting pseudo-observables, e.g. rH->.4f, analytic continuation is performed 
only after integration over all variables but the Higgs virtuality. The role of pseudo-observables is then very 
similar to what one had for Lep physics, e.g. the hadronic and leptonic peak cross-sections 

<^!l = 12.^, af^l2.I|^, (50) 

where the POs are fully inclusive. Already at that time it was clear that the expression for realistic observables 
(RO) at arbitrary s requires a careful examination because of non triviality in off-shell gauge invariance, see 
Rcf. I60i. 
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To summarize, our proposal aims to support the use of Eq. (j27p for producing accurate predictions for the 
Higgs hneshape in the SM; Ea. (l?7l) can be adapted easily for a large class of extensions of the SM. Similarly 
Eg. ([5^ should be used to define the Higgs-boson production cross-section instead of ctos or ctofs of Eq. ([551) . 

To repeat an obvious argument the zero-width approximation for the Higgs boson is usually reported in 
comparing experimental studies and theoretical predictions. The approximation has very low quality in the 
high-mass region where the Higgs boson has a non negligible width. An integration over some distribution 
is a more accurate estimate of the signal cross-section and we claim that this distribution should be given 
by the complex propagator and not by some ad hoc Breit-Wigner. 

As we have already mentioned most of this paper is devoted to studying the Higgs-boson lineshape in the 
high-mass region. However, nothing in the formalism is peculiar to that application and the formalism itself 
forms the basis for extracting pseudo-observables from experimental data. This is a timely contribution to 
the relevant literature. 



3.2 Production cross-section 

Before presenting detailed numerical results and comparisons we give the complete definition of the 
production cross-section; let us define (^ = zs, k = vs, and write 



prod = ^ I PDF af;tn = E V ("X"' 



all(C,'«,m,/^F), (51) 

. . ~ — .1 , u 

where zq is a lower bound on the invariant mass of the H decay products, the luminosity is defined by 

f ^ dx / V \ 

Ctj{v) ^ — fi{x,fiF) fj [-,fJ-F) , (52) 

where fi is a parton distribution function and 

1 Ck r*j°t(C) 

Therefore, tTy_j.H+x(C5 Mr) the cross section for two partons of invariant mass k (z < u < 1) to produce 
a final state containing a H of virtuality = z s plus jets (X); it is made of several terms (see Ref. |5lj for a 
definition of Act), 



y2 <^i3~^H+x{(, K, /iR, Mf) = 0-gg_>H (5 ( 1 - - ) + - (ACTgg^Hg + ACTqg^Hq + ACTcjq_^Hg + NNLO) . (54) 



As a technical remark the complete phase-space integral for the process pi +pj — > pfc + {/} {pi — XiPi etc.) 
is written as 

J d^.j^f = J rf^prod J d$dcc = J d^pkS+{pl) Yl d^<li6^i<lf)S'^ i^Px+P] -Pk -^q}j 

= fd'kd^QS+ipl)S'{p,+p,~p,-Q) I W d^qi5+{qT)5^[Q-Y.qi\ (55) 



where / rf$dcc is the phase-space for the process Q — > {/} and 

rf$prod = S [ dz J d^Pkd^Q 5+{pI) 5 (Q2 _ Q 0(Q„) ^4 (^^ + _ _ g) 

dzdvdi f d^pkd^Q 5+{pI) 5 (Q2 _ Q) e{Qo) 5^ {p, + ~pk-Q) 
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(56) 



Egs. dSip and (1531) follow after folding with PDFs of argument and Xj^ after using xi = x, Xj — v/x and 
after integration over i. At NNLO there is an additional parton in the final state and five invariants are 
need to describe the partonic process, plus the H virtuality. However, one should remember that at NNLO 
use is made of the effective theory approximation where the Higgs-gluon interaction is described by a local 
operator. Our variables z,v are related to POWHEG parametrization [5S], F,^, hy Y = ln{x/y/z) and 



As already mentioned perturbation theory breaks down in for the complex pole when = 1-74 TeV 
(around 930 GeV for on-shell mass). When we consider the lineshape for = 900 GeV (Fos — 647 GeV) 
in a window of, at least, ±2 Fos around the peak we go above 2 TeV in the Higgs virtuality; it is worth 
noting that for /in = 1-5(2) TeV the on-shell width is 3.38(15.8) TeV; these results follows from using 
Prophecy4f [¥7] which includes two-loop leading corrections and it is worth observing that the LO result 
at /iH = 2 TeV is 3.7 TeV with a 58% one-loop corrections and an estimated 270% two-loop contribu- 
tion/uncertainty. 

These considerations give additional support to the GPP - scheme as the only viable option for defining 
the Higgs signal in the very high mass region; indeed the F^*(s) in Eq. ipT]) is replaced by a Fh*(sh) in the 
CPP - scheme. The OFFP-scheme remains the most pragmatic alternative but the tail of an heavy Higgs 
boson lineshape extending above 1.3—1.5 TeV should be taken with due caution. 

4 The issue of gauge invariance 

As discussed in Section [2] only the CPP-scheme of Eq . (|46|) should be used and the background should be 
included in the calculation. Let us consider the pragmatic OFFP-scheme of Eq. (|45p in more detail: here we 
use the Higgs propagator with its complex pole but production and decay are computed at arbitrary Higgs 
virtuality and not at the complex pole. As far as LO production is concerned, e.g. the ggH one-loop fermion 
triangle, there is never an issue of gauge-parameter dependence in going off-shell; in this respect higher order 
QCD corrections are not a problem. 

Consider now the decay, i.e. F^' in Eg. ([5^ : the amplitude v4(H — > F), for each final state |F > and 
as long as we include the complete set of diagrams at one-loop order, is gauge-parameter independent if the 
Higgs boson is on its mass-shell. However, as soon as we put an external leg off-shell, the amplitude must be 
coupled to the corresponding physical source I, i.e. A{1 — >■ H — ^ F), and only the complete process I — ^ F is 
gauge-parameter independent. The latter does not exclude the existence of subsets of diagrams that satisfy 
the requirement but this can only be examined case-by-case. To rephrase it, if the Higgs boson is off shell, in 
LO and NLO QCD in most cases the matrix element still respects gauge invariance, but in NLO EW gauge 
invariance is lost, unless the CPP-scheme is used. 

How to deal with the OFFP-scheme? Technically speaking, we have a matrix element 



where s is the virtuality of the external Higgs boson, /in is the mass of internal Higgs lines and Higgs 
wave- function renormalization has been included. The following happens: /(sh , sh) is gauge-parameter 
independent (CPP-scheme) to all orders while /(/^h , /^h) i^ gauge-parameter independent at one-loop but 
not beyond, /(s , ^^) is not. In order to account for the off-shellness of the Higgs boson we defined the 
OFFP-scheme by choosing (at one loop level) /(s , s), i.e. we intuitively replace the on-shell decay of the 
Higgs boson of mass /j,h with the on-shell decay of an Higgs boson of mass -^/s and not with the off-shell 
decay of an Higgs boson of mass ^h- The same applies for the NLO EW correction to production. 
The proof that CPP-scheme satisfies gauge-parameter independence can be sketched as follows: 



v = z/{l-0- 




(57) 



S{s) 



^prod(sH) ^dcc(sH) 




(sh)[1 + 5{jhGsh)] 



-1/2 



= 0, 



(58) 



[1 + 5{jh(^h)] (s-sh)' 



-proi 
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where ^ is an arbitrary gauge parameter, ^'^^(sh) is the derivative of ShhIs) computed at s = sh and 
^prod,dcc is the amphtude for production/decay including vertices and wave-function renormahzation con- 
stants for the external lines. Note that Eg . ([55)) follows from the use of Nielsen identities, see Ref. [BTj. An 
example is discussed in Appendix El 

At one loop the renormalized Higgs self-energy, computed in the i?^ -gauge, can be written as 

sill - , = 4h - ,^^l) + {s- Ml) ASlll (,s , Mh , e) • (59) 

Finally, to clarify the statement that /(/i^ , /i^) may violate gauge invariance beyond one loop we consider 
the H — )■ YY decay. The LO amplitudes depends on the bare Higgs mass and standard renormahzation, 
introducing the on-shell mass through the real part of the one-loop self-energy, leads to a violation of the 
Ward-Slavnov- Taylor identities (WSTI), see Refs. [62 63 64 ; this will happen at two- loop above the WW 
threshold. Therefore, the real part of the Higgs self-energy stemming from mass renormahzation must be 
traded for the complex expression, even if the external Higgs boson is assumed to be an on-shell particle. 
Eq. (|58p gives an additional argument against Breit-Wigner distributions: only the propagator 

' Zhh = 1+ ^""^^^'^""^^"\ (60) 



reproduces the correct factors to be combined with production and decay vertices. 

To conclude the discussion we summarize few well-known points which, however, are not always taken into 
proper account. Consider the partial decay width H — WW, hiding the fact that also the Ws are unstable 
objects; one way of constructing it is to consider the amplitude for the process and derive at tree-level the 
following expression: 



A 



H^WW 



Next we compute the imaginary part of the one-loop H self-energy due to internal charged lines. We obtain: 



At(s) = -l^(l-'-^)(l-H!^r^ (62) 

Unitarity of the theory requires proportionality between F^^^^ and Im . This is never a problem for 
the fermion final states but we see already the paradox in the bosonic sector: from one side we compute 
H —J' WW, from the other all charged bosonic lines (physical or not) circulate in the loop. It is trivially true 
that 



= r^'^^''^^(Mi) = Im SUMll (63) 



a result that satisfies unitarity and gauge invariance; however, going off-shell is different. There are matters 
of principle that do not matter in practice: if we consider a 'light' Higgs, say below 200 GeV, the difference 
between the on-shell width computed in the -gauge and 7h is tiny and, aside from the neighborhoods of 
unphysical thresholds {£_ = M^/(4M-^),^ — A/y/(4M|)), negligible over a wide range of values of the gauge 
parameter ^ |31) . However, consider a typical implementation of the H production cross-section: 



(.-M2)2+M2Fh(M2)- 



Z n r2T^ /..^2^ ' (64) 



13 



What it is meant by rH(M^) in the denominator of Eq. (j64p ? Since we are deahng with the denominator 
of the propagation function for a Higgs boson of virtuahty s and since we would like to respect unitarity, it 
has to be interpreted as the imaginary part of the self-energy. At one loop, for s = gauge invariance 
is not violated. But what about the rH(s) that appears in the numerator? Again, because of unitarity, 
it should be interpreted as the imaginary part of the off-shell self-energy, i.e. the self-energy for a H of 
mass Mfj and virtuality s; however, this quantity has a f -dependent part which is a function of the Higgs 
virtuality. Therefore, no matter how small is the on-shell width for a light Higgs boson, there will be an 
s and f -dependent distortion of the lineshape. Of course there will never be a problem for a calculation 
assembling the whole process in one stroke; however, defining the signal is another story. 

We conclude this Section by showing explicitly the delicate point in extracting the pseudo-observable 
decay width from the complete amplitude. We have 



Sis) 



1 



■ShhIs) - 5'hh(sh) 



^prod(s) 



1 



^prod(s) 



S - SH 

5'hh(s) - "ShhIsh) 



S - SH 



S - SH 
-1/2 1 



l^dcc(sH) + ^doc(s) - ^dcc(sH) I 



^dcc(s) 



S - SH 



•ShhIs) - S'hh(sh) 



S - SH 



-1/2 



= <5'prod(s) 



S - SH 



Sdcc{s), 



(65) 



where the matrix element for the decay can be decomposed into 



5dcc(s)= l + ixHH(s) 



Xhh{s) = S'hh(sh) 



■ShhIs) - 5'hh(sh) 



SH 



— 1/2 ^ 
-^H ^dcc(sH) + 1 - iHH(s) Adoc(s) - Adcc(sH) 



rHH(s) 



S'hh(s) - S'hh(sh) 



(66) 



The decay width is defined in terms of Z^^^^ ^doc(sH) while the rest is O (s — sh) 



5 Numerical results 

In the following we will present numerical results obtained with the program HTO (G. Passarino, unpub- 
lished) that allows for the study of the Higgs-boson- lineshape, in gluon-gluon fusion (ggF), using complex 
poles. HTO is a FORTRAN 95 program that contains a translation of the subroutine HIGGSNNLO written 
by M. Grazzini for computing the total (on-shell) cross-section for Higgs-boson production (in ggF) at NLO 
and NNLO |65I66I67] (HTO has a library for one-loop integrals that extends the one of Ref. [55]). It is 
worth noting that the ggF production mechanism dominates up to /in = 700 GeV but for higher masses 
the vector-boson fusion (VBF) mechanism starts do compete, with a ratio ggF/VBF of 1.11 at 1 TeV. The 
following acronyms will be used: 

FW Breit-Wigner Fixed Width (Eq.(l35])) 

RW Breit-Wigner Running Width (Eq. (p4)) ) 

OS parameters in the On-Shell scheme 

Bar parameters in Bar-scheme (Eq. p6|) ') 

FS QCD renormalization (factorization) scales fixed 

RS QCD renormalization (factorization) scales running 
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Table 3: The re- weighting factor w for the (total) integrated Higgs lineshape, (jP^p/ct'^^ using 7h ~ in Eq. (|32p . 



MH [GeV] 


[GeV] 


[pb] 


[pb] 


crP™P [pb] 


w 


200 


1.43 


5.249 


5.674 


5.459 


0.962 


300 


8.43 


2.418 


2.724 


2.585 


0.949 


400 


29.2 


2.035 


1.998 


1.927 


0.964 


500 


68.0 


0.8497 


0.8108 


0.7827 


0.965 


600 


123.0 


0.3275 


0.3231 


0.3073 


0.951 



All results in this paper (but those in Table [HI) refer io ^/s = 7 TeV and are based on the MSTW2008 PDF 
sets [69]. For complex W, Z poles we use Eq.((9|) with HXSWG standard input. For the t complex pole we 
use [TOj 



OS 



with a LO on-shell width 



8V27r 



M3 



(67) 
(68) 



with an on-shell mass Aft = M^'^ = 172.5 GeV. 



The integrand in Ea. (j5ip has several peaks; the most evident is in the propagator and in HTO a change 
of variable is performed, 



c = 



1 



1 + im/my 



+ 7h tan (/^hTh z') 



{am + au) C- ar, 



arctan , cm = arctan 

MH7H 



Similarly, another change is performed 



— exp{(l — u^) In z}, X = exp{(l — y^) Inv}. 



(69) 



(70) 



In comparing the OFFBW-scheme with the OFFP-scheme one should realize that there are two sources of 
difference, the functional form of the distributions and the different numerical values of the parameters in 
the distributions. To understand the impact of the functional form we have performed a comparison where 
(unrealistically) 7h = is used in the Higgs propagator; results are shown in Table [3] 

In Table [D we compare the on-shell production cross-section as given in Ref. [1. with the off-shell cross- 
section in the OFFBW-scheme (Eq.((32])) and in the OFFP-scheme (Eq.((l5|)). The total Higgs boson width, 
present in Eg. ([251) . is computed by interpolating the tables of Ref. [T]. 

An important question regarding the numerical impact of a calculation where the on-shell Higgs boson 
is left off-shell and convoluted with some distribution is how much of the effect survives the inclusion of 
theoretical uncertainties. In the OFFBW-scheme we can compare with the results of Ref. [21], at least 
for values of the Higgs boson mass where their results are available, i.e. below 300 GeV. We compare 
with Tables 2—3 of Ref. [Tl at 300 GeV and obtain the results shown in Table [5] where only the QCD scale 
uncertainty is included. At these values of the Higgs boson mass both on-shell production and off-shell 
production, sampled over a Breit-Wigner, are compatible within the errors. The OFFP-scheme gives a 
slightly higher central value of 2.86 pb. The comparison shows drastically different results for high values of 
the mass where the difference in central values is much bigger than the theoretical uncertainty. 

In Table [6] we show a more detailed comparison of the production cross-section in the OFFBW-scheme 
between iHixs of Ref. [211 and our calculation. A is the percentage error due to QCD scales and 6 is the 
percentage ratio HTO/iHixs. 
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Table 4: Comparison of the on-shell production cross-section as given in Ref. jTj with the ofT-shell cross-section in 
the OFFBW-scheme {Eq.^) and in the OFFP-scheme (Eq. (|45)l '). 



MH[GeV] 




7H[GeV] 




a«W[pb] 


crP™P[pb] 


500 


68.0 


60.2 


0.8497 


0.8239 


0.9367 


550 


93.1 


82.8 


0.5259 


0.5161 


0.5912 


600 


123 


109 


0.3275 


0.3287 


0.3784 


650 


158 


139 


0.2064 


0.2154 


0.2482 


700 


199 


174 


0.1320 


0.1456 


0.1677 


750 


248 


205 


0.0859 


0.1013 


0.1171 


800 


304 


245 


0.0567 


0.0733 


0.0850 


850 


371 


277 


0.0379 


0.0545 


0.0643 


900 


449 


331 


0.0256 


0.0417 


0.0509 



Table 5: The production cross-section in pb at /in = 300 GeV. Result from HTO is computed with running QCD 
scales. 

Tab. 2 of Ref. [J Tab. 3 of Ref. [I] Tab. 5 of Ref. [21] HTO RS-option 

9^9+0.14 f, .r+0.16 9 r7+0.15 9 oi +0.25 
^-^^-O.IS ^■^'^-0.22 ' -0.22 ^■°-^-0.23 



Table 6: Comparison of the production cross-section in the OFFBW-scheme between iHixs (Table 5 of Ref. [21]) and 
our calculation. A is the percentage error due to QCD scales and 5 is the percentage ratio HTO/lHlxs. 





CTiHixs [pb] 


AiHixs [%] 


o-HTO [pb] 


Ahto 


[%] 


5[%] 


200 


5.57 


+7.19 -9.06 


5.63 


+9.12 - 


9.30 


1.08 


220 


4.54 


-1-6.92 -8.99 


4.63 


+8.93 - 


-8.85 


1.98 


240 


3.80 


-1-6.68 -8.91 


3.91 


+8.76 - 


8.51 


2.89 


260 


3.25 


-1-6.44 -8.84 


3.37 


+8.61 - 


-8.22 


3.69 


280 


2.85 


+6.18 -8.74 


2.97 


+8.49 - 


7.98 


4.21 


300 


2.57 


+5.89 -8.58 


2.69 


+8.36 - 


7.75 


4.67 
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Table 7: Scaling factor OFFP/OFFBW schemes (Eqs.(|3Sl) and ([331) for the invariant mass windows of Ea. (f72)) . 





n -4 


n = —3 


n = -2 


n = -1 


n = 


n = 1 


n = 2 


n = 3 


500 


1.24 


1.13 


1.08 


1.23 


1.32 


1.11 


0.97 


0.88 


600 


1.39 


1.20 


1.08 


1.27 


1.40 


1.12 


0.94 


0.84 


650 


1.50 


1.24 


1.08 


1.30 


1.44 


1.13 


0.93 


0.81 


700 


1.62 


1.28 


1.09 


1.34 


1.49 


1.13 


0.92 


0.79 


750 


1.78 


1.34 


1.10 


1.41 


1.57 


1.14 


0.90 


0.77 


800 


1.99 


1.41 


1.12 


1.49 


1.64 


1.14 


0.89 


0.75 



Another quantity that is useful in describing the lineshape is obtained by introducing the peak mass, 

Mp^oak=(A^g')'+m', (71) 

an by considering the following windows in the invariant mass, 

Mpeak + ^ , A/peak + rgS] , n-0,±l,.... (72) 

In Table[7]we present the ratio OFFP/OFFBW for the invariant mass distribution in the windows of Eg. ([7^ . 

In Table [5] we include uncertainties in the comparison; both OFFBW and OFFP are computed with the 
RS option, i.e. running QCD scales instead of a fixed one. Note that for a{iJ,) we define a central value 
(Jc — o^iJi) and a scale error as [a^ , a^] where 



cr = mm 



MS 



[m/2,2m] CT(Ai), a+ = max^g[jj/2,2jr] ^(Ai), (73) 



and where fJ. — ^J,n — /ip and Jl is the reference scale, static or dynamic. 

In Table ini we compare results from Table 5 of Ref. [H] with our OFFBW results with two options: 1) 
fip are fixed, 2) they run with C/4. For three values of reported we find differences of 2.0%, 3.7% and 

4.7%, compatible with the scale uncertainty. Furthermore, we use Ref. T for input parameters which differ 

slightly from the one used in Ref. [21]. Note that also the functional form ot the Breit-Wigner is different 

since their default value is not the one in Eq.(|34|) but 



BW(.) = i y^^P^ , (74) 



TT 



S — 



where rH(-\/s) id the decay width of a Higgs boson at rest with mass i/s. 

In Table [TU] we use the OFFP-scheme of Eg. ipS)) and look for the effect of running QCD scales in the 
H invariant mass distribution. Differences are of the order of 2—3% apart from the high-mass side of the 
distribution for very high H masses. As we have mentioned already, the only scheme respecting gauge 
invariance that allows us for a proper definition of pseudo-observables is the CPP-scheme of Ea. (|15)) . It 
requires analytical continuation of the Feynman integrals into the second Riemann sheet. Once again, 
the definition of POs is conventional but should put in one-to-one correspondence well-defined theoretical 
predictions with derived experimental data. 

In Table [TT] we give a simple example by considering the process gg — >■ H at lowest order and compare 
the traditional on-shell production cross-section, see the l.h.s. of Eg. ([55)1 . with the production cross-section 
as defined in the CPP-scheme. Therefore we only consider gg H aX LO (i.e. v ^ z) and put H on its real 
mass-shell or on the complex one. Below 300 GeV there is no visible difference, at 300 GeV the two results 
start to differ with an increasing gap up to 600 GeV after which there is a plateau of about 25% for higher 
values of ^h- 
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Table 8: Production cross section with errors due to QCD scale variation and PDF uncertainty. First entry is the 
on-shell cross-section of Table 2 of Ref. p], second entry is OFFBW (Eq.JiaJ), last entry is OFFP Ea. (|45]l . 



MH[GeV] 


CT[pb] 


Scale [' 


%] 


PDF [%] 


600 


0.336 


+ 6.1 - 


5.2 


+ 6.2 


- 5.3 




0.329 


+ 11.2 - 


11.3 


+ 4.7 


- 4.8 




0.378 


+ 9.7 - 


10.0 


+ 5.0 


- 3.9 


650 


0.212 


+ 6.2 - 


5.2 


+ 6.5 


- 5.5 




0.215 


+ 12.4 - 


12.2 


+ 5.1 


- 5.3 




0.248 


+ 10.2 - 


10.0 


+ 5.3 


- 4.2 


700 


0.136 


+ 6.3 - 


5.3 


+ 6.9 


- 5.8 




0.146 


+ 13.9 - 


13.3 


+ 5.6 


- 5.8 




0.168 


+ 11.0 - 


11.1 


+ 5.5 


- 4.6 


750 


0.0889 


+ 6.4 - 


5.4 


+ 7.2 


- 6.1 




0.101 


+ 15.8 - 


14.5 


+ 6.1 


- 6.3 




0.117 


+ 12.2 - 


11.9 


+ 5.8 


- 5.1 


800 


0.0588 


+ 6.5 - 


5.4 


+ 7.6 


- 6.3 




0.0733 


+ 18.0 - 


16.0 


+ 6.7 


- 6.9 




0.0850 


+ 13.7 - 


13.0 


+ 6.1 


- 5.5 


850 


0.0394 


+ 6.5 - 


5.5 


+ 8.0 


- 6.6 




0.0545 


+ 20.4 - 


17.6 


+ 7.2 


- 7.5 




0.0643 


+ 15.7 - 


14.3 


+ 6.5 


- 6.0 


900 


0.0267 


+ 6.7 - 


5.6 


+ 8.3 


- 6.9 




0.0417 


+ 22.8 - 


19.1 


+ 7.7 


- 8.0 




0.0509 


+ 18.1 - 


16.0 


+ 7.0 


- 6.6 



Table 9: Comparison of results from Table 5 of Ref. with our OFFBW results (Ea. (f43ll 'l with two options: 1) 
/in , /ip are fixed, 2) they run with (^/4. 





iHixs 




HTO 










MH[GeV] 


f7[pb] 


Scale [%] 


CTi [pb] 


Scale [%] 


(72 [pb] 


Scale [%] 




220 


4.54 


+ 6.9 - 9.0 


4.63 


+ 8.9 - 8.9 


4.74 


+ 7.3-8 


.7 


260 


3.25 


+ 6.4 - 8.8 


3.37 


+ 8.6 - 8.2 


3.49 


+ 7.9-8 


.6 


300 


2.57 


+ 5.9 - 8.6 


2.69 


+ 8.4 - 7.8 


2.81 


+ 8.4 - 8 


.5 



Table 10: Scaling factor in the OFFP-scheme (Ea. (|45|) ') running/fixed QCD scales for the invariant mass distribution. 
Here \x,y\ = [A/p.^j, + a: ^g^ Mpeak + V Tg^] • 



AiH[GeV] 


[-1,-5] 




[-i,o] 




U ' 2J 


[5.1] 


600 


1.023 


1.021 


1.019 


1.018 


1.017 


1.016 


700 


1.027 


1.023 


1.021 


1.019 


1.018 


1.019 


800 


1.030 


1.024 


1.021 


1.021 


1.021 


1.136 
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Table 11: 7h and Fg as a function of fin- We also report the ratio between the LO cross-sections for gg — )■ H in the 
CPP-scheme and in the OS-scheme. 



Mh[GcV] 7h[GcV: 



r^=^[GcV] acPp(sH)/aos(A/^i'") 



300 7.58 

350 14.93 

400 26.66 

450 41.18 

500 58.85 

550 79.80 

600 104.16 

650 131.98 

700 163.26 

750 197.96 

800 235.93 

850 277.00 

900 320.96 

950 367.57 

1000 416.57 



8.43 0.999 

15.20 1.086 

29.20 1.155 

46.90 1.183 

68.00 1.204 

93.10 1.221 

123.00 1.236 

158.00 1.248 

199.00 1.256 

248.00 1.262 

304.00 1.264 

371.00 1.263 

449.00 1.258 

540.00 1.252 

647.00 1.242 



The production cross-section in Ea. (|54p is made of one term corresponding to gg ^ _ff which we denote 
by a'^ and terms with at least one additional jet, cr^ . The invariant mass of the two partons in the initial 
state is K = V s with z < w < 1, C being the Higgs virtuality. In Table [T^ we introduce a cut such that 



and study the dependence of the result on Vc by considering a* — a'-* + and 6 = cr*/cr*. The results show 
that a cut z < v < l.Ql z gives already 85% of the total answer. 

LHC is now running with a center of mass energy of 8 TeV; in Table [T^ we have shown a comparison for 
the production cross-section at 7 TeV and 8 TeV. 

On the left-hand side of Figure [1] we compare the production cross-section as computed with the OFFP- 
scheme of Eq.^ or with the OFFBW-scheme of Eq.l^ for = 600 GeV. For the OFFBW-scheme 
we use Breit-Wigner parameters in the OS-scheme (red curve) and in the Bar-scheme of Ea.(p6| (blue 
curve). Deviations from the OFFP-scheme are maximal in the OS-scheme and much less pronounced in the 
Bar-scheme. 

On the right-hand side of Figure [T] we show the effect of using dynamical QCD scales (invariant mass of 
the di-photon system) for the gg-^H + X^-X + yy cross-section at /in = 400 GeV. 

On the left-hand side of Figure [2] we consider the on-shell production cross-section (t'-*^(pp H) (black 
curve) which includes convolution with PDFs. The blue curve gives the off-shell production cross-section 
sampled over the (complex) Higgs propagator while the red curves is sampled over a Breit-Wigner distribu- 
tion. The observed effect is substantial even in the low mass region. 

On the right-hand side of Figure [2] we show differential if factors (ctnlo/clo etc.) for the process pp — ?> 
(H — >• 4e) -I- X, comparing the fixed QCD scale option, /ir — fip — Mh/2, and the running scale option, 
Mr = fJ-F — Af (4e)/2. For running QCD scales the i^T factor is practically constant over a wide range of the 
Higgs virtuality. 

On the left-hand side of Figure[3]we show the normalized invariant mass distribution in the OFFP-scheme 
with running QCD scales for Higgs-boson masses of 600 GeV (black), 700 GeV (blue), 800 GeV (red) in the 
windows Mpcak ± 2ros where Mpcak is defined in Eq. (j7ip . 

On the right-hand side of Figure |3] we show the normalized invariant mass distribution in the OFFP- 
scheme (blue) and OFFBW-scheme (red) with running QCD scales for a Higgs-boson mass of 800 GeV in 
the window Mpeak ± 2 Fos- 



z < V < min{wc z , 1} 



(75) 
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Table 12: The efTect of introducing a cut on the invariant mass of the initial-state partons. 



Mh[GcV] 




a"[pb] 


(T-' [pb] 


cr*[pb] 


6[%] 


600 


no cut 


0.2677 


0.1107 


0.3784 






1.001 




0.0127 


0.2804 


74 




1.01 




0.0473 


0.3149 


83 




1.05 




0.0918 


0.3595 


95 


700 


no cut 


0.1208 


0.0469 


0.1677 






1.001 




0.0054 


0.1262 


75 




1.01 




0.0201 


0.1409 


84 




1.05 




0.0391 


0.1599 


95 


800 


no cut 


0.0632 


0.0218 


0.0850 






1.001 




0.0025 


0.0657 


77 




1.01 




0.0094 


0.0726 


85 




1.05 




0.0182 


0.0814 


96 



Table 13: Comparison of the production cross-sections in the OFFP-scheme (Eq.(|45j) at 7 TeV and 8 TeV. 



AiH[GeV] 


cr [pb] 7 TeV 


cr[pb] 8 TeV 


ratio 


600 


0.378 


0.587 


1.55 


650 


0.248 


0.392 


1.58 


700 


0.168 


0.271 


1.67 


750 


0.117 


0.194 


1.67 


800 


0.0850 


0.145 


1.77 


850 


0.0643 


0.113 


1.77 


900 


0.0509 


0.0922 


1.87 



Finally, in Figure|4]we show the normalized invariant mass distribution in the OFFP-scheme with running 
QCD scales at 600-800 GeV in the window Mpeak±2ros for 7 TeV (blue) and 8 TeV (red). The invariant 
mass distribution in the OFFP-scheme with running QCD scales for a Higgs-boson mass of 800 GeV (in the 
window Mpeak ± 2ros) is shown in Figure [SJ the blue line refers to a center-of-mass energy of 8 TeV, the 
red one to 7 TeV. 

5.1 Results in CPP-scheme 

Here, we will briefly summarize what was covered in Section [3. II for the CPP-scheme. referring to Eq. ()53p 
we can say the following: the production amplitude, including QCD corrections, is always ^-independent 
and we can keep full off-shellness; the decay amplitude, beyond LO, is not and should be evaluated at the 
Higgs complex pole in the CPP-scheme. To give an example of the results we select the process gg — 
H — >■ Z'^Z'^, where the two Z -bosons in the final state are also taken at their complex pole, thus giving a 
full gauge-invariant quantity. We will show results in the OFFP-scheme (Eq. (j45|) ) and in the CPP-scheme 
(Eq. p^ ). Caveat emptor: one should not interpret their difference as a source of theoretical uncertainty 
since comparisons should only take place for the full S ® B matrix element and what we are doing in the two 
schemes simply amounts to shifting some part from the known S-amplitude to the unknown S/B-interference. 
Furthermore, only the CPP-scheme is a fully consistent recipe to define a gauge-invariant signal; the only 
reason why we presented most of our results in the OFFP-scheme is that it is much simpler to implement in 
any preexisting calculation while the CPP-scheme requires analytical continuation into the second Riemann 
sheet. 

In Figure [8] we show the invariant mass distribution in the OFFP-scheme (black) and in the CPP-scheme 
(red) for = 700 GeV (left) and = 800 GeV (right) in the window [Afpoak-Mpeak + 2ros] for the 
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Table 14: Comparison of the production cross-sections at 8 TeV: OFFBW-scheme, OFFP-scheme, on-shell results 
of Ref . [71] and the Breit- Wigner approach of Ref. [2^ . 





f7[pb] OFFBW 


cr[pb] OFFP 


cr[pb] Ref. [71, 


cr[pb] Ref. [22] 


300 


3.74 


3.87 


3.32 


3.45 


280 


4.10 


4.25 


3.67 


3.82 


260 


4.60 


4.79 


4.16 


4.34 


240 


5.29 


5.53 


4.82 


5.05 


220 


6.21 


6.43 


5.71 


6.00 


200 


7.48 


7.53 


6.95 


7.32 


190 


8.31 


8.35 


7.78 


8.21 


180 


9.42 


9.45 


8.85 


9.42 


170 


10.59 


10.62 


10.11 


10.72 


160 


12.24 


12.24 


11.76 


12.66 


150 


13.88 


13.88 


13.57 


14.43 


140 


15.78 


15.79 


15.59 


16.53 


130 


18.05 


18.05 


18.04 


19.14 


125 


19.41 


19.41 


19.49 


20.69 



process gg — >■ H — >■ Z'^Z'^. The shoulder in the low mass side of the resonance (in the CPP-scheme) is due 
to the fact that the H -width is very large, the production cross-section is increasing for decreasing invariant 
mass and there is no suppression from the decay which is taken at the Higgs complex pole. Thus, the 
decrease in the distribution for small invariant masses is delayed by the large width. The high invariant 
mass behavior reflects again a constant (at the complex pole) decay amplitude. 

5.2 Results in the low Higgs mass region 

For lower values of the Higgs mass we compare the OFFBW-scheme and the OFFP one (based on Eq.((T7])) 
at 8 TeV in Table [TU here we also show the on-shell results of Ref. [7T] and the Breit- Wigner (running width 
scheme) approach of Ref. The results of Table [Til show a good convergence of OFFBW/OFFP results 
to the on-shell ones for low values of the Higgs mass. 



6 QCD scale error 

The conventional theoretical uncertainty associated with QCD scale variation is defined in Eg. ([75]) . Let 
us consider the production cross-section at 800 GeV in the OFFP-scheme with running or fixed QCD scales, 
we obtain 

RS 0.08497 FS 0.07185+^-^^°. (76) 

The uncertainty with running scales is about twice the one where the scales are kept fixed. One might wonder 
which is the major source and we have done the following; in estimating the conventional uncertainty in the 
RS-option we keep fj,p—JI^ — (/4: and vary between Ji/2 and 2/1. The corresponding result is 

RS 0.08497+^;^% only. (77) 

Therefore, the main source of uncertainty comes from varying the factorization scale. However, in the 
fixed-scale option we obtain 

FS 0.07185+^;^^°, m only (78) 

and /ip -variation is less dominant. In any case the question of which variation is dominating depends on the 
mass-range; indeed, in the ONBW-scheme with fixed scales at 140 GeV we have 

ONBW 12.27 ^10.1%, m and /ip 
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ONBW 12.27 till: m only (79) 

showing /iR dominance. 

In order to compare this conventional definition of uncertainty with the work of Ref. |72j we recall their 
definition of p% credible interval; given a series 

k 



a = ^c„<, (80) 



n—l 

the p% credible interval a ± Acr^ is defined as 

AaP = max ({c}) p%, p% < 

Uc nc + 1 

r -1 — 1/nc rj 

AcjP -a^+imax({c}) (n, + l)(l-p%) , p%> (81) 

L J -\~ 1 

with Uc = k — I + 1. Using this definition we find that the 68% (90%) credible interval for /iR uncertainty 
in the production cross-section at 800 GeV is 3.1%(5.3%) in substantial agreement with the conventional 
method. For ONBW-scheme at 140 GeV the 90% credible interval for uncertainty is 8.0%. It could be 
interesting to extend the method of Ref. |72] to cover fip uncertainty. 



7 Residual theoretical uncertainty 

As we mentioned before the accuracy at which the imaginary part of sh can be computed is not of the 
same quality as the next-to-leading-order (NLO) accuracy of the on-shell width. It would we desirable to 
include two- and three-loop contributions as well in 7h and for some of these contributions only on-shell 
results have been computed so far. Here we follow the analysis of Ref. [TH] where the authors try to improve 
the solution of the complex pole equation by using the available on-shell information. It is again a method 
based on some expansion and does not include the full SM, QCD corrections and complex poles for internal 
W, Z lines; it amounts to consider only the Higgs-Goldstone Lagrangian of the SM. Nevertheless it is very 
useful to give a rough estimate of the missing orders. The authors of Ref. [18] define a double expansion of 
the self-energy, 

Suh{s)=AMI+B{s~M^) + j^{s-M^)' + ... A = ^a„(^^)", (82) 
etc. In this way we can write 

Gf ( Gf Ml a„ \ 

7H - — ^ ai 1 + — ^ ^ + . . . . (83) 

The ratio 02/01 has been computed and we can estimate that the first correction to 7h is roughly given by 

Gp /i?r 

fe= 0.350119 ^ . (84) 

Changes in 7h range from 2.3% at 400 GeV to 9.4% at 750 GeV. In general, from Ref. [18 we do not see 
very large variations up to ITeV with a breakdown of the perturbative expansion around when 1.74 TeV. 
Let us elaborate a bit more: in the Higgs-Goldstone model one has 

— = 1.1781 5H + 0.4125 5^ + 1.1445 <?3 g^ = ^L^ (85) 
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Table 15: Theoretical uncertainty on the production cross-section, the height of the maximum, the position of the 
half-maxima and the area of the resonance. 



/iH[GeV] 


Sh[%] 


^crPiod 


[%] 




AC-, AC+[GcV] 




SA[%] 


600 


5.3 


-6.0 H 


-6.3 


-10.8 +11A 


(-2.5, -f2.5) (-1-2.5, - 


-2.5) 


-4.8 +6.0 


700 


7.2 


-8.0 H 


-8.6 


-14.8 +16.0 


(-9.0, -f8.0) (+4.0, - 


-4.0) 


-7.0 +11.8 


800 


9.4 


-9.7 + 


10.6 


-19.3 -f21.5 


(-18.2, -fl8.2) (+6.1, 


-6.1) 


-8.7 +9.5 



giving 7H = 168.84(LO), 180.94(NLO), 186.59(NNLO) GeV for /in = 700 GeV. Using the three known terms 
in the series we estimate a 68% credible interval of 7h = 186.59 ± 1.93 GeV. The difference NNLO — LO is 
17.8 GeFand in the full SM our estimate is 7h = 163.26+ 11.75 GeV. Therefore, using 7h (1 ± (5h) we can 
give a rough but reasonable estimate of the remaining uncertainty on the lineshape. Of course, one could 
plot the lineshape with an error band (an example is shown in Figure |6] for /xh = 700 GeV) but it is better 
to quantify the uncertainty at the level of those quantities that characterize the resonance. Referring to 
Eq. ([5T|) we introduce the following quantities: 

(jP™'^ total production cross-section 

^^prod 

I] = differential distribution 

{Cmax , Smax} thc maximum of the lineshape 

{C± , — Smax} thc half-maxima of the lineshape 

A = / d(T, the area of the resonance between half-maxima (86) 

where C is the Higgs virtuality. Changing 7h to 1 ± S^i gives a variation on the production cross-section 
which is linear in i5h, as shown in Table [T51 similarly, the variation in the peak height goes as 2^h with no 
appreciable variation in its position, Cmax- In Table [15] we also show the variation (in GeV) in the position 
of the half-maxima and in the area of the resonance. 

To summarize our estimate of the theoretical uncertainty associated to the signal: 



• the remaining uncertainty on the production cross-section is typically well reproduced by ((5h + 1 ) [7o 

• 5],„ax changes approximately with the naive expectation, 2 Jjj [%] 

• The shift induced in the position of the half-maxima are more complicated but a rough approximation 
is given by 

AC±«±^^7Hfe, M±=/iJj/2 (/iH±7H)'/'. (87) 

How to use these results depends on the specific analysis, H WW Ivlv, Ivqq and H ^ ZZ ^ lllL 
etc. For example, the lineshape information is used in H — > WW — Ivqq analysis [73]. The H — ZZ — )• llvv 
analysis use transverse mass spectrum for the llvv system. If the analysis is a simple "counting experiment" , 
therefore not taking the Higgs lineshape into account, the uncertainty on the area and the half-maxima can 
be neglected. Therefore, informations on the shape of the resonance will become more relevant once more 
advanced analysis techniques (based on differential distributions as discriminant variable) will be used. 

The factor ry*(C) in Eg. ip?)) deserves a separate discussion. As explained in Section S] it represents the 
"on-shell" decay of an Higgs boson of mass y/( and we have to quantify the corresponding uncertainty. The 
staring point is r*°* computed by Prophecy4f [4^ which includes two-loop leading corrections in GpM^, 
where Mh is now the on-shell mass. Next we consider the on-shell width in the Higgs- Goldstone model. 
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Table 16: Theoretical uncertainty on the total decay width, F^* in Ea. (|27p . Fp is the total width computed by 
Prophecy4f, AF gives the credible intervals and 5F[%] is the ratio Ar/Gammap. 





rp[GeV] 


(5r[68%] 


5r[95%] 


600 


123 


0.25 


0.42 


700 


199 


0.62 


1.03 


800 


304 


1.35 


2.24 


900 


449 


2.63 


4.38 


1000 


647 


4.72 


7.85 


1200 


1205 


13.1 


21.7 


1500 


3380 


34.7 


57.8 


2000 


15800 


98.9 


165 



discussed in Refs. |18|74j . We have 

^L = E-V' = ^»- A.^. ,88, 

Let Fp — Xp ^ the width computed by Prophecy4f, we redefine the total width as 

F ^ 

n=0 

where now oq = Xp — Xhg- As long as A is not too large we can define a p% < 80% credible interval as 
(following from 02.3 < ai) 

rtot(C)=rp(C)±AF, AF = ^max{|ao |,ai}p%A4v/C- (90) 

It is easily seen that for ^/C — 929 GeV the two-loop corrections are of the same size of the one-loop 
corrections and for >/C = 2.6 TeV one-loop and Born become of the same size. Not that the expansion 
parameter in Eq.((88|) is one for = 1.57 TeV. In Table [TBI we present both the 68% and the 95% credible 
intervals. This should be compared with the corresponding estimate in Ref. |2] of 17%(My/l TeV) for 
Mh > 500 GeV. Due to the large effect induced by radiative corrections for large values of Mh and to 
the fact that the two-loop correction is only known in the Higgs-Goldstone Lagrangian we have devoted 
Appendix [B] to a detailed discussion of the validity of the whole approach. 

It is clear that it does not make much sense to have an error estimate beyond 1.3 TeV and, therefore, 
all results for the Higgs lineshape that have a sizable fraction of events in this high-mass region should 
not be taken too seriously. Here, once again, the only viable alternative to define the Higgs signal is the 
CPP-scheme. 

To summarize: above 0.93 TeV perturbation theory becomes questionable since the two-loop corrections 
start be become larger than the one- loop ones; above 1.3 TeV the error estimate also becomes questionable 
since the expansion parameter is A = 0.7 and the 95% credible interval (after inclusion of the leading two- 
loop effects) is 32.2%. The invariant mass distribution for /in = 800 GeV and the corresponding uncertainty 
introduced by Ftot(C) ^-I's shown in Figure [71 

The theoretical uncertainties coming from 7h and Ftot can be combined but a meaningful prediction re- 
quires that we cut the Higgs virtuality {() at some upper value for which we have selected Cmax = (1-5 TeV)^. 
Results are shown in Table [171 and in Table [TSl otherwise the non-perturbative increasing of Ftot(C) over- 
compensates the decreasing of aij^^i with a negligible smearing from the invariant mass distribution. For 
completeness we present the leading (in the limit — > 00) if -factor for the H — > VV amplitude (V — W, Z) 
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Table 17: Total theoretical uncertainty on the production cross-section at 7 TeV. The total is obtained by considering 
the THU on 7h and on Ttot with a cut ^ < 1.5 TeV. 



/iH[GeV] 




'^[%] 


600 


-5.5 


+5.9 


700 


-7.0 


+7.5 


800 


-7.7 


+8.8 


900 


-7.0 


+8.9 



Table 18: Total theoretical uncertainty on the production cross-section at 8 TeV. The total is obtained by considering 
the THU on 7h and on Ttot with a cut < 1.5 TeV. 







[%] 


600 


-5.2 -i 


-5.8 


700 


-6.3 -i 


-7.4 


800 


-6.7 -f 


-8.4 


900 


-6.0 -f 


-8.1 



in the OS-scheme and in the CPP-scheme: 



GfM?, / GfM^ . 2 
OS if - 1 + ai ^ + a2 ^ ^ 



16V27r2 V 16v/27r2 - 
where the coefficients are [73] 

ai = 1.40- 11.35i 02 = -34.41 - 21.00 i (92) 
Above 1 TeV the NNLO term dominates the if -factor. 



7.1 Interference signal/background 

The most important factor that has to be observed here is that our analysis of the residual theoretical 
uncertainty refers to the signal only. In the current experimental analysis there are additional sources of 
uncertainty, e.g. background and Higgs interference effects. As a matter of fact, this interference is partly 
available and should not be included as a theoretical uncertainty; for a discussion and results we refer to 
Refs. [49150] . In particular, from Refs. [15] we see that at /in = 600 GeV (the highest value reported) the 
effect is about +40% in the window C = 440—560 GeV , is practically zero at the peak and reaches —50% 
after C, — 680 GeV (no cuts applied). For the total cross-section in gg — )• Ivl'v' at /xh = 600 GeV the effect 
of including the interference is already +34% and rapidly increasing with /in . 

We stress that setting limits without including the effects of the interference induces large variations in 
rate and shape that will propagate through to all distributions. Therefore, any attempt to analyze kinematic 
distributions which arc far from the SM shape may result in misleading limits. 



8 Conclusions 



In the last two decades many theoretical studies lead to improved estimates of the SM Higgs boson total 
cross section at hadron colliders [75l76l51l77l78l79l80t5T] (including electroweak NLO corrections [82183] ) 
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and of Higgs boson partial decay widths [84147) . Less work has been devoted to studying the Higgs boson 
invariant mass distribution (Higgs-boson-hneshape) |21I22I55] . 

In this work we made an attempt to resolve the problem by comparing different theoretical inputs to 
the off-shellness of the Higgs boson. There is no question at all that the zero-width approximation should 
be avoided, especially in the high-mass region where the on-shell width becomes of the same order as the 
on-shell mass, or higher. 

The propagation of the off-shell H is usually parametrized in terms of some Breit-Wigner distribution; 
several variants are reported in the literature, e.g. Eq. (p4| and Eq.(|74|). see Ref. [57] for a discussion. We 
have shown evidence that only the Dyson-resummed propagator should be used, leading to the introduction 
of the H complex pole, a gauge- invariant property of the 5 -matrix. 

Finally, when one accepts the idea that the Higgs boson is an off-shell intermediate state the question of 
the most appropriate choice of QCD scales arises. We have shown the effect of including dynamical choices 
for QCD scales, /xr,/zf, instead of a static choice. 

For many purposes the well-known and convenient machinery of the on-shell approach can be employed 
but one should become aware of its limitations and potential pitfalls. Although the basic comparison between 
different schemes was presented in this work, at this time we are repeating that differences should not be 
taken as the source of additional theoretical uncertainty; in Section |4] we have shown that the CPP-scheme, 
described in Section l3T] (see Eg. (j46|) ) . is the only the only theoretically consistent tool for describing the 
lineshape; it requires, however, the computation of Feynmann diagrams on the second Riemann- sheet. As a 
consequence, our simple pragmatic recommendation is to use the OFFP-scheme described in Section [3TT] (see 
Eq. pS)) ) with running QCD scales. The associated theoretical uncertainty has been discussed in Section [7] 
where we also examine the limitations induced by an apparent breakdown of the perturbative expansion in 
the ultra- heavy Higgs boson limit. 

Finally, the most severe problem faced by a critical appraisal of the heavy Higgs boson lineshape is the 
missing inclusion of interference with the background, especially in the light of gauge invariance issues. 
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Figure 1: In the left figure we show a comparison of production cross-section as computed with the 
OFFP-scheme of Eq.(|45|| or with the OFFBW-scheme of Ea. (|43)l . The red curve gives Breit-Wigner 
parameters in the OS-scheme and the blue one in the Bar-scheme of Eq.((36)|. In the right figure we show 
the effect of using dynamical QCD scales for the production cross-section of Eas. (|51ll - (l53p . 




Figure 2: In the left figure the blue curve gives the off-shell production cross-section sampled over the 
(complex) Higgs propagator while the red curves is sampled over a Breit-Wigner distribution The black 
curve gives the on-shell production cross-section. The right figure shows differential if factor for the 
process pp — >■ (H — >■ 4e) + X, comparing the fixed scale option, /ir — fip = Myi/2, and the running scale 
option, /iR = fiF = Af(4e)/2. 
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Figure 3: The normalized invariant mass distribution in the OFFP-scheme with running QCD scales 
(left) for 600 GeV (black), 700 GeV (blue), 800 GeV (red) in the windows Mpoak±2ros- The normahzed 
invariant mass distribution in the OFFP-scheme (blue) and OFFBW-scheme (red) with running QCD 
scales (right) at 800 GeV in the window Afpoak i 2ros- 
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Figure 4: The normalized invariant mass distribution in the OFFP-scheme with running QCD scales for 
600 GeV (left), 800 GeV (right) in the windows Mp^ak ± aFos- The blue hue refers to 8 TeV, the red 
one to 7 TeV. 
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Figure 5: The invariant mass distribution in the OFFP-scheme with running QCD scales for 800 GeV 
in the window Afpg^j; ± 2 Fos- The blue line refers to 8 TeV, the red one to 7 TeV. 
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Figure 6: The invariant mass distribution in the OFFP-scheme with running QCD scales for fiu = 
700 GeV in the window Mpeak i Fos- The red lines give the associated theoretical uncertainty as 
discussed is Section [T] 
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Figure 7: The invariant mass distribution in the OFFP-scheme for /in ~ 800 GeV in the window 
[Mpgai;— Mpcak + 2ros] with the error band due to the theoretical uncertainty introduced by r^*((^) 
discussed is Section [T] 
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200 400 000 700 SOO 1000 1200 200 000 800 1000 ir>00 

Figure 8: The invariant mass distribution in the OFFP-scheme (black) and in the CPP-scheme (red) for 
= 700 GeV (left) and fin = 800 GeV (right) in the window [A'/p^ak - 2ros - Mp^ak + ZFos] for the 
process gg — ^ H — )> Z^Z". 
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Figure 9: The resonant part of the process ij — > H + fc where i,j and k are partons (g or q). We distinguish 
between production, propagation and decay, indicating the explicit dependence on Mandelstam invariants. 
If s denotes the pp invariant mass we distinguish between vs, the ij invariant mass, zs, the Higgs boson 
virtuality and t which is related to the Higgs boson transverse momentum. Within NR we include the 
Non-Resonant contributions as well as their interference with the signal. NR' includes those terms that 
are of higher order in the Laurent expansion of the signal around the complex pole. At NNLO QCD the 
parton A; is a pair k\—k2 and more invariants are needed to characterize the production mechanism. 
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A Appendix: Nielsen identities at work 

In this Appendix we consider one specific example of application of Nielsen identities ( see Ref. [5T]) at 
one loop level, namely we explain of to get gauge-parameter independence in the process 

R{P)^Zip,,^i) + Zip2,ly). (93) 

We work in the i?j -gauge; for any quantity /(^) we write a decomposition 

/(e)-/(l) + A/(0, A/(1) = 0. (94) 

The most important quantity is the Higgs self-energy, 

SuH{s)^slll + 0{g*)^^EMs) + 0{g'), s = -P\ (95) 
Let A/jj be the renormalized Higgs mass, we obtain that the one-loop, ^ -dependent, part can be given as 

As(j^^(e, s, M^) = (s - Ml) (e, s, M^). (96) 
The main equation is the one for the Higgs complex pole, 

SH - m2 + sill ' , Ml) = 0, (97) 

from wich we derive 

m2=,h+0(5'), ^^'^(f, sh, .sh) -0. (98) 
Next we consider the one-loop vertices contributing to H — >■ ZZ and obtain a one-loop amplitude 

4'^(H^ZZ)= (^i^)<5^, + i;(i)p2MPi.)e^(pi,Ai) e''(p2,A2). (99) 
Using the decomposition 

we obtain the following results: 

1. after reduction to scalar form-factors there are no scalar vertices remaining in AvJ^'' , sh , sh)- 

2. Furthermore, AV^'^'' (C, sh , sh) 0. 

Both results are expected since the corresponding terms could not cancel against anything else. Furthermore 
we compute renormalization Z -factors for the external legs [53] 

/-c . „ ^ c(i) 



^-Ml + sill , ^ , Ml) = {s - ,Sh) [l + ^™ (^-^-^H)-^HH(e.^H,.H) - 

= (i + Zh) (s-sh) +0((s-sh)') . (101) 



It follows that 

AV^^ (C, SH , sn) - \l AZh(0 + AZ^iO] = 0, (102) 



where the lowest-order is defined by 



^^°^ = -5^, -o4 = ^, Ml^, = s^,s,. (103) 
At one loop it is enough to evaluate all functions at s = but only sjj is ^ -independent. 
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B Appendix: equivalence theorem for (off-shell) virtual bosons 

In this Appendix we discuss the extension of the equivalence theorem to virtual vector-bosons. Formally, 
the theorem states, at the 5 -matrix level, that 

S[V^]={tC)'"S[4>], (104) 

where n is the number of external, longitudinal, vector-bosons and C = 1 at tree level. In Ref. [85 it is 
shown how one can get C = 1 at higher orders through a choice of the </> renormalization constant. We 
consider the amplitude for the process 

H(P) ^ f (gi) + f (fci) + f fe) + f (fc2), (105) 

where all fermion are massless. The amplitude can be split according to the number of resonant vector- 
bosons, 

A^Aav + Aiv + Aov, (106) 
where A2Y starts from LO in perturbation theory while Aiy,Aqy start from NLO. Furthermore we have 

A,y = [F,S^'^ + F,Pj^] A^Jpi)A^^(p2) J"(gi,fci) J^(<?2,fc2), (107) 

where Ay is the vector-boson propagator and J"'^ are (conserved) fermion currents. In the ^ = 1 gauge 
there are no scalar propagators involved and each V-propagator is proportional to Sf^a etc. 

We are interested in the leading term of the amplitude in the limit Affj — > 00. This term derives from 
a balance between H -couplings (0 stands for a generic Goldstone-boson) and H -propagators; at one- loop 
(or higher) Aiy is given by boxes and Aiy by pentagons. Since external fermions are massless we can easily 
see that there are not enough internal scalars to contribute to the leading term. 

We will use the formalism of Refs. |86I87I88] and will write 

S,.-^^ei^Mel,Ap)+ E e^^ip,X)elJp,X). (108) 

^ A=-l, + l 

Since the fermion currents are conserved we can replace the numerator in the V-propagators with the sum 
over polarizations. A convenient choice for the polarizations is the following: 

eL^,{pl) ^ -Nl (pi ■ P2Pif, + siP2f,), eLp(p2) = -A^l (Pi ' P2m + 52 Pi > (109) 
where A^i,2 are the normalizations, p^ — —Si and 

e±M,X) = ^ [n^{p^)+^XN^{p^)], N^{p^) = (s.)"'/' e^^a/jp ""(ft) e^^(K)P^ (HO) 

nt^iPi) = iN]_ef,ai3pk" p^p^, nf,{p2) = i Nlsf^aiSp p^ p^i- (111) 
With this choice one obtains 

E e^Mb'^)el^(p,A) = (5^^ - - eL^(p)e£^(p). (112) 

A=-l, + l ^ 

The first question to answer is about L -polarization dominance in Higgs decay for Mh — cxj^j. For the Fd, Fp 
parts in Ea. (ll07p we have to evaluate 



E e{pi Ai) • e{p2 X2) e* (pi, Ai) e^(p2, A2), 



^We gratefully acknowledge A. Denner and S. Dittmaier for an important discussion on this point. 
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^ e(piAi) -^2 6(^2 A2) •pie*(pi,Ai)e^(p2,A2). 



(113) 



Ai,A2 



First we introduce invariants for the complete process 

qi- ki 



-si, q2-k2 = --S2, qi-q2 = -^S3, 



qi-k2 = -^S4, ki-q2 = -^S5, ki-k2 = -^se, 



(114) 



where the invariants satisfy s, = M^. Therefore, not all scales in the problem can be small as compared 
to M^. For the V squared-propagators we have 



(S - M2 )2 + r% M2 My Fv 



(S-M2)2 



+ ... 



(115) 



Therefore, always in the limit Mh — >■ 00, we can assume that si^2 ~ My 0. In the rest-frame of the 
Higgs boson and using massless fermions in the final state, the condition si_2 =0 requires Sj = M^/A for 
i = 3, . . . , 6. Next we introduce 



D 



L;L 



eL(pi) • eL(p2) , -D_L_;_L_ = e_L(pi,-l) • e_L(p2,-l) 



etc. in the limit Si,2 — >■ and Mr — >■ 00. We obtain 

1 I S1+S2 



D 



L;L 



Mi| + 0(1), 



4 si,S2 2 .S1.S2 

while all other combinations are finite in the limit. We also introduce non-diagonal elements, 

Nxapr = e{pi, A) • e{p2, a) [e{pi, p) ■ e{p2, r)] *, 



(116) 



(117) 



(118) 



where e{p, 0) = ei,{p) and e{p, ±1) = e±{p, ±1). Also these elements are finite or zero in the limit si, S2 — ^ 0. 
Furthermore, we introduce other quantities 



Pt,. 



L:L 



eh{pi) ■P2 ehip2) -Pi 



etc. (including non-diagonal terms) and find 



16 S1S2 4 S1S2 



(119) 



(120) 



while all other combinations are finite in the limit. Finally we define the fermion part of the amplitude, 
squared and summed over spins: 



rL= ^ |j(g,fc)-eL(p) , ri;A = ^ U(9,fc)-ei(p,A) 



where p = q + k and 



We obtain 



J''(9, k) = u{k) Y {V+ 7+ + V- 7-) v{q), 7± = 1 ± 7'- 



(121) 

(122) 
(123) 
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where s — —p^. We also compute non-diagonal elements 

TA.a = '^('?' ^) • ^) [-^(^^ fc) • e(p, (t)] \ (124) 

spins 

and find 

To^±i^2V2i{V^-Vl) s + 0{s^) , r^i^+i^2{Vf + V^) s + 0{s^) . (125) 

Therefore, all F factors are suppressed in the low-s region and, according to Eqs. (lll7p - (jl20l) . only longitu- 
dinal vector bosons contribute to the leading term for A/h — oo. In conclusion, for the double-resonant part 
of the amplitude, only the longitudinal polarizations matter. One-loop corrections are computed exactly 
in Ref. |47] which also includes two-loop leading corrections; they are extracted from the Higgs-Goldstone 
model and Prophecy4f multiplies the complete LO width with the NNLO correction factor. For the con- 
tributions of the transverse vector-bosons it is not correct but irrelevant in the limit Mjj oo, as we have 
explicitly shown. 

Before concluding our argument we need to explain how the equivalence theorem (see Ref. [SS]) works 
in practice for our off-shell vector bosons. The theorem is based on a Ward-Slavnov- Taylor identity which 
relates the amplitude for a V -boson contracted with its momentum to the corresponding amplitude where 
the V -boson is replaced by a Goldstone-boson 0. However, in our case we have 

eL/.(Pl) = -^L {Pl •P2PlA< + SlP2p^), (126) 

and the longitudinal polarization is not simply proportional to the momentum. To see how it works in the 
Mh —J' oo limit, we consider the LO HW^W~ vertex contracted with longitudinal polarizations 

Vl;L = -5^% eL(?3i) • eL(p2), bl p.ipi) = -NI {Lpi ■ p2Pi^_, + R sip2 ^_,) , (127) 

where we have artificially introduced the L, R coefficients. It is easy to show that only the part of Vl : l 
survives in the limit, giving 

Vl-^l-'^9^MI (128) 

2 ^/SlS2 

which, after identification of Si with Af^ is exactly minus the II(/)+(/<^ vertex, as dictated by the equivalence 
theorem (but only in the limit A/h ^ oo). As expected, the WSTI relating 9^V^ to Myip, is valid only for 
on-shell particles. The same holds for the combination 

- g Mw ehipi) ■ P2eLip2) ■ Pl, (129) 

showing that, in the limit — ^ oo, we can extract the leading behavior from the Higgs-Goldstone model 
and we can multiply the LO width with the NNLO correction factor, channel-by-channel. The only delicate 
point concerns the WSTI which substitutes two (contracted) V -lines with (/) lines. When only one of the two 
is contracted and the other is multiplied by the wrong momentum there is an additional term in the WSTI. 
We illustrate this singly-contracted WSTI at LO in Figure [TUl if the J -source is physical (pj ■ Jfj.{P2) = 0) 
then the sum of the first two diagrams is zero, otherwise there is a remainder given by diagram c) of Figure llOl 
where Faddeev-Popov ghosts are involved. A similar situation happens also at higher orders, as shown in 
Figure [TTl where the dash- arrow line represents Faddeev-Popov ghosts and the off-shell V -source contracted 
with the wrong momentum can absorb a V -ghost pair. One can argue, on the basis of the explicit calculation, 
that this term is never leading at one-loop but we have not found a convincing argument that this is the 
case to all orders. 

The above arguments, strictly speaking, apply to the decay of an on-shell Higgs boson. If this decay 
is part of a complete process, e.g. gg — ?> 4f (where we have Higgs- resonant diagrams) then the following 
happens: let g be the SU{2) coupling constant and A = g'^ M^/M^; the work of Ref. [85] shows that, for 
g^ /X — > 0, the equivalence theorem is true to lowest nonzero order in g an to all orders in A. However, there 
is no reason to expect the theorem to hold when using a propagator with a Breit-Wigner width, which is 
equivalent to summing a subset of diagrams of higher order in A. 

We conclude the Appendix with another comment: technically speaking Ea. (|104p requires that S[(j)\ is 
computed in the full SM, not only in the Higgs-Goldstone model. We have verified in simple examples that 
the difference is subleading when A/h — >■ oo. 
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J''{P2) r{P2) r{p-2) 



Figure 10: Example ol a LO simply-contracted WSTI where the dash-hne represents the Higgs boson, 
solid line represent vector-bosons and dash- arrow-lines represent Goldstone-bosons or Faddeev-Popov 
ghosts. If the source J is physical then a + b — c — 0, otherwise a + b + c — 



Pi 




XP'I 



Figure 11: Example of a remainder for a NLO simply-contracted WSTI where the dash-line represents 
the Higgs boson, sohd line represent vector-bosons and dash-arrow-hnes represent Faddeev-Popov ghosts. 
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C Appendix: complex poles, real kinematics and pseudo-observables 



In this Appendix we discuss how to implement complex poles in the amphtude for H — > F within the 
CPP-scheme. For sake of simphcity we consider the process H(P) — > Z(pi) + Z(p2) where both Zs are on-shcU 
while the Higgs boson is off-shell, i.e. = — s. At LO the -matrix element is 

Mz 

Q^ipi) ■ e{p2). (130) 

cos t) 

Squaring the amplitude and summing over polarization introduces a which is a kinematic factor and 
should be replaced with s 

2 

^ I^LO =V2Gfs^1-Axz + 12xI), xz = — . (131) 

spins 

At LO there is never a problem with gauge parameter dependence. At NLO there is an amphtude 

^NLO ' ^^h) = s Frf (s , ) + (s , Ml) p'i pl^ , (132) 

where the arises from vertices and propagators, always with P^ — —s. Renormalization in the complex- 
mass scheme requires the replacement — 7> sh ; taking the residue at the complex pole also requires s — > sh 
and guarantees gauge parameter independence. When we go to the 5 matrix element and compute the NLO 
correction we obtain 

2Re^ ^io^NLO = -144(y2GF)'^'Mws'[xz^^^^^ + -^(l-i— + ^xz) ReFp], (133) 

S D Xz 'J 

spins 

where the form factors are computed at s = M^ — sh- When there are more particles in the final state the 
analytical continuation of the form factors should be done according to Eqs.(P5 )) - P^ and generalizations, 
followed by the integration over the phase space with real angles and energies. 

C.l Analytic continuation 

We continue this Appendix emphasizing, yet once more, that complex kinematics is not needed since we 
need to perform analytic continuation in the S* -matrix element (because of gauge invariance) while keeping 
the phase-space real; this is the approach introduced in Ref. |89| . the so-called leading-pole-approximation 
(LPA). In Appendix A of Ref. [52] it is shown that, in order to have a pole on the second Riemann sheet 
independently of the other momentum variables, analytic continuation of the scattering amplitude in the 
variable is required (fc being the momentum of the unstable particle). Therefore the 5 -matrix has to be 
analytically continued in around the positive energy cut to the second sheet and the state vector for the 
corresponding unstable particle must not depend on /c*^; QFT provides an explicit construction for this state 
as described in details in Ref. [52] . 

Consider the process where an off-shell Higgs boson of virtuality decays into n massless particles of 
momentum g^; then the quantity of interest is 



rH^„(C) = / rf$i^„ 

H^n {{qi ■ q,}) 



1 
1 



2 



s.c 

2 



where {/} is the set of 3n — 7 invariants, linearly independent, needed to describe the process. Another 
quantity of interest is the production cross-section for an off-shell Higgs boson plus jets, e.g. g{pi) +g(p2) ^• 
g('7i) +H((72), with ql = C, 



2 
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^'Isn-l Ehgg-H+„g({/})f, (135) 



An extension to complex kinematics is possible and the relevant work can be found in Ref. j52l; here we 
repeat the construction for completeness. We consider distorsions of the original real contour that project 
one-to-one onto the original contour under the projection that sends complex momentum onto its real part. 
The contour of integration may be distorted in any bounded region provided the integrand is analytic through 
the region of distorsion. 

Suppose that, in a given process, kc is the momentum of an unstable )incoming or outgoing) Higgs 
particle and k is its projection back to the real axis; the general form of kc will be (here we work in the 
Bjorken-Drell metric) 



^h-^A.h(1 + :^) -;^7h(i + ^) . (136) 
The vectors p, q must satisfy 

-q^ = mL P-q = - 2 MH7H- (137) 







kl = 


Mh 








1/2 


1 






71 


" Mh ] 





To have a resonance k must be a timelike momentum and we can go to a frame k — yjQ , j where ( is the 
Higgs virtuality. We have chosen the following analytic continuation, yielding the complex pole: 

Stepl) ki^^ L{kQ,k) k, fcc==(VC,o), (138) 

where L is a Lorentz boost; 



Step 2) k^ km = , Oj , (139) 
where the analytic continuations is defined by 

fcA = ^((2-A)C + AAi|,0,0,2A(M^-C)), 0<A<1. (140) 

Finally, we boost back in the k direction and continue into the complex plane as follows: 

Step 3) -^kc=p + iq, p = (rc k° , k) , q=(/mfc^,o). (141) 



The continuation in Eas. (|138l) - (|139p is the real, on-shell, projection in LPA. Since q^ > from Eq. (|137l) it 
follows that p^ > 0. Any such state may be transformed by a real Lorentz transformation into one with 
standard momentum k defined by 



( P'l n n lG{p,q) 



where G(j), q) — p^q^ — {p ■ q)^ is the Gram determinant of p, q. 

Once complex boosts are allowed states with momentum kc satisfying fc^ — fi^ can be constructed from 
Eq. (|142l) . However, states appearing in scattering processes have a timelike momentum. Therefore, the only 
unstable particle states that appear in physical processes have a momentum kc that becomes timelike when 
continued back to the real axis, i.e. > is needed to produce resonances while the associated q is not 
constrained. In the additional two cases we have standard momenta (i.e. in each case kc is related by a real 
Lorentz transformation to the standard momentum k 

2 . - / AV + (p.g)^ AV-(p-g)^ A _ 
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9 < 0, p 



g2 



0, 



(143) 



The possibility of having two parametrizations connected by a complex boost (an clement of L(C) the 
complex, homogeneous Lorentz group) may seem like posing a question of uniqueness, but it does not since 
scalar products are defined without complex conjugation of one of the complex four-vectors. 

Since we perform our analytic continuation on the invariants (see Eqs. (ll34|) - (|135|) ) we have to define the 
mapping. Consider again the process g(pi) +g{p2) g(fci) + H(fc2) and define s = {pi +^2)^, t = {pi — fc2)^ 
and 



t 



1 - 



( 1 — cos ( 



t) 



) (1 - cos( 



(144) 



where s has been kept real to preserve s — 4i?^, {E being the c.m.s. energy) and 9 is the scattering angle, 
also real. We add one condition to the mapping {/} — >■ {Ic}, that after continuation all the invariants are 
analytic functions of sh ■ It is easily seen that the condition is fulfilled by 



p + iq, p= 2^ (s + ^l, 0, 0, -s + /i|^ , g = - (1,0,0,1), 



2^s 



i.e. = 0, corresponding to 



(sH + S , , , SH - s) ■ 



(145) 



(146) 



In principle there is no need to discuss crossed channels separately. Indeed, any function A{{p}), after going 
through complex momenta, may be analytically continued back to the real axis but with negative energies 
for some of the incoming or outgoing particles. The result of this analytical continuation is the amplitude 
for a different physical process (e.g. H — >■ ggg). One should only pay attention to the fact that 



A{{I})=<a\S\b> A{{I})^<b\S\a> (R}) ^ 4 ({/J) , 



(147) 



for complex invariants. Let us consider the decay of a Higgs boson of virtuality = ( into thre massless 
particles of momenta ki . . . k^. We define invariants Sij = {ki + kj)"^ and angles Oij. Select one invariant, say 
S12, to remain real and keep ^12 also real. The analytic continuation is defined by 



'23 



C-S12- ((C-Si2)'-4cCsi2 



1/2 



'13 



+ ((C- Sl2)^ -4CCS12) 



1/2- 



(148) 



with c = cot^(0i2/2). Of course, we can keep any of the real and define analytic continuations of the same 
function, rH^3, in different portions of the 3 -dimensional, complex, space Si2,si3 and S23: in the sense of 
analytic continuation there is no ambiguity. The choice of which invariant remains real is equivalent to the 
choice of which momentum is fixed by conservation. To follow the argument we write 



r 



H->3 



(C) 



{{I}) 



2VC 



rfsi3 ^ L4h^3 (C , S12 , 513^ 



-7= / dsi2 / d COS 6*12 J (cos 6*12 , si2 , C) y^ K^H-^3 (C: Sl2 008 6*12; 

'VC Jo J-1 I 



(149) 



where 9 is the condition (C — 512)^ > 4cCsi2, 



1 /■+^ I 
rH^3(C) = ^ / ^0086*12 / ds 12 J {cos 9 12 , s 12 , C) y2\A}i ^3 {C, s 12 cos 9 12) 
z vC J-i Jo „ , I 



^ , , (712 1712 

C- = ( cot — - CSC — 



C- 



(150) 
(151) 
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A more democratic approach in dealing with H decay into n massless states is defined by the foUowing choice: 



Pi^ Pi{\, k) ^ [X + in) pi-^ p1 = {a + i(i) Pi, z = l...n + l (152) 
where pf = Q for i > 1 and 

a^-r = — , a/3 = -- — —. 153 
Pi 2 p{ 

This continuation respects momentum conservation and mass-sheh condition for (massless) stable particles 
and works as long as the continuation passes through analytic regions of the integrand as A varies from 1 to 
a and k, varies from to /3. 

There is never any problem with fully inclusive quantities but a careful discussion is needed if we want 
to consider differential distributions with cuts: in this case one would like to integrate a function, whose 
arguments have been analytically continued into the complex plane, over a real (fixed) phase-space. In this 
case rH^3(C) becomes 

—= / dcose'i2 / dsi2 J(cos6'i2 , si2 , C) y~lr^H^3 (SH , Sl2 cos6'i2) (154) 
^VQ J-1 Jo , J 

To close the discussion we should say that the definition of signal is always conventional since differences 
(that are non-resonant) will always be included into background plus interference. Furthermore, when main 
emphasis is on differential distributions in pp — > 4f there is little meaning in splitting the whole process 
into components; vice versa, when the problem is with extracting pseudo-observables, e.g. rH^4f , analytic 
continuation is performed only after integration over all variables but the Higgs virtuality. Nevertheless, if 
one wants to introduce cuts on differential distributions we have two alternatives. We can define!^ 

rH^3(sH) =Rer^^3(sH), (155) 

1 f~^^ f ^ I 2 

rH^3('SH) ^ ry , / dcOs6'i2 / dsi2 J(cOs6li2 , Sl2 , C) 5Ir^H^3 (SH , Si2 C0S6'12) (156) 



where 7 is a curve in the complex S12 -plane connecting the origin with the pont S12 = Sjj (sjj being the 
continuation of ^_ from (. to sjj) provided the integrand is anlaytic through 7. Finally, we project back into 
the real S12 -axis with a change of variable 

( dsi2= [ (isi2 >/c (si2 , Ci sh) , (157) 
Jj Jo 

where Jc is the Jacobian of the transformation. 

Alternatively the NLO corrections to a process with n + 1 external legs can be written as a sum over D, 
the set of scalar one-loop, n -I- 1 -leg diagrams 

^|AH^„({/})|'=Re ^ Qd(C,{/}) / [dxrf] l^,-^--^ ({x4 , C, {/}) , (158) 

where Qd are functions of the invariants, {x} are the Feynman parameters that parametrize the loop inte- 
gral for diagram d, Sd is a simplex in the a; -space, Vd is a quadratic (linear) form in the Feynman parame- 
ters (invariants) and is a positive integer (4 — e is the space-time dimensionality). We can introduce cuts 
on the real kinematics, 

(159) 



rr>„(C) - ^ / E |^H^„ ({/}) 



We gratefully acknowledge S. Pozzorini for an important suggestion on this point. 
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perform first the integration over the phase space with cuts, perform analytic continuation in one variable 
(C ~^ ■5h) while doing the integration over Feynman parameters (which requires continuation into the second 
Riemann sheet). The functions Qd may contain square roots, depending on the masses of the final states 
and on their number. For massless final states and n < 3 each Qd is a rational function of the invariants. 
For higher values of n it is well-known that phase-space integrals may lead to elliptic functions. 

To give an example we consider H(p) — g(/ci) -|- g(fc2) -|- g{kj,) with s = [pi — ki)^ and t = {p— fca)^. For 
instance, there will be a contribution from the direct box with internal top quarks; the fully extrapolated 
contribution corresponding to Q = 1 from Eq. (|158|) can be cast into the following (manifestly finite) form: 



dx / dy 
Jo 



Li2 I ^° ) + Li2 (r-a) - Li2 (n) 



1 



{1- x)y {x - y) 



ra = ^v{l-x), n^-^y{l-y). 

If we impose a cut, e.g. (1 — zi)C^ < s < {1 — Z2) C, the result is equally compact, 
f-i 



(160) 
(161) 



dx I dy 
^0 



Li2 yzi 
- Li2 ( Z2 



n - Tg 

l-ra 

n - Tg 

1 - Tg 



Li2 yzi 

Li2 I Z2 



I - ZlTg 

Tg 

I - Z2rg 



- Li2 yzi 

+ Li2 I Z2 



n 



I - ZlTg 

n 

I - Z2rg 



{l~x)y{x-y) 



(162) 



The technique we have described is easily extendible to an arbitrary number of legs. Consider a process 
^ with all incoming momenta; a convenient way of describing the boundaries of the phase space, 
the relations between vectors and invariants and the non-linear constraints that arise when > 6 id the 
following. If pi is the momentum associated to the Higgs boson we redefine n-dimensional momenta and 
introduce components 

Pi = (m 7 0„_i) , (163) 



where 0^ is a string of k zero components. Then we introduce 



hi = 1, 



(164) 



'^kl 

"'kl ; rs 

etc. The vector p2 is defined by 



l^kl 



P2 



Next n-dimensional vectors will be 

0„-3: 



P3 



Pi 



hi. 



with row k and column I removed, 
with row r and column s removed. 



M2 



, 0„ 



det/i2 
det/ii 



det 



'Mi 



0„-4 . 



Vdet det 
det h% . 32 




Vdet /ii det ' Vdet det ' V det 



det 



(165) 
(166) 



(167) 



etc. If A'^ = 5 then momentum conservation gives P5 ^ — Pi- If, instead TV = 6, pg follows from momentum 
conservation whereas p^ requires a fifth component; if n = 4 the vanishing of the fifth component requires 
det/i^ = 0. The boundary of the phase space is determined by the equations 



det /i^ < 0, 



, det h 



N-l 



< 0. 



(168) 
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As we have mentioned several times the analytic continuation requires knowledge of the branch cuts of the 
integrand. When the Higgs boson is the only massive external line a box can be expressed in terms of 
i? -functions (see Appendix B of Ref. 90 ): 



R{ya , 2/1) = / dy 



y-yo 



ln(y-2/i)-ln(yo-yi) 



(169) 



where yo,i depend on the invariants of the problem. The function R has the following branch cuts in the 
complex yo, yi -space: 

1) yi = (1 - x)yo 

2) yi = {^-x)yo + x 

3) 2/0 = 2; yi 

4) yo = xyi + 1- X 
Pentagons etc. are always linear combination of boxes. 



< a; < 1 



C.2 Pseudo-observables 

We consider now the problem of extracting pseudo-observables from matrix elements. First we introduce 
the two-body phase-space, 

$2(s,si,S2) = ^^ / n d^P.S+{pl + s,) 5\P-p,~p2) 

1 Al/2(s,Si,S2) 



16 TT-y/s 



s-(yir+v^)^ 



16 TT 



with = —s. Consider again the process H(P) — !■ Z(pi) + Z(p2) with an amplitude given by 



4° 



w 



cos^ 



,41 - g 

Afw 



(170) 



(171) 



The form factors are function of the external invariants and of the internal masses but they are computed 
at the Higgs complex pole, i.e. 

(172) 



so that they are gauge parameter independent. We introduce the following quantity: 
D2 = Y^^^^2 (s,M|,Af|) \ {A% + Al,)e^{pi)el{p2)\\ 



(173) 



which is the product of a Higgs propagator times the S' -matrix element for the Higgs boson into a pair of 
longitudinal Z bosons. We obtain 



Ml 



Th^z^z, (s) 
f2 nj2 



VV2 (.,Af|,Af| 



9 



1 Xis,Mi,Mi) 



1 S - Ml 



cos^ 6 4 A'/| cos^ I 



cos^ 9 



\{s,Mi,Mi)Fp^ — 



4 A/4 -v-'-z,-z;-P 2 



ReTi), 



\{s,MlMl) 



A/| 



(174) 



(175) 



There is no problem with gauge parameter dependence in Ea. (ll75p : nevertheless, it is useful to separate the 
residue of the Higgs complex pole and also to introduce Z complex poles obtaining a pseudo-observable 



16 7r^| SH 



Re r, 







^ Ti), |sh|»|sz 



cos^ 9 



(176) 
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where T is obtained with the substitution s — > sh and Af| — > sz in T. As a second step we consider 
H ^ ZZ e+e~|i+n^ and define 



Da = 



spins 



1 



Az(si) PI Az(s2) P 



where J is the (conserved) leptonic current and 



is the phase-space for H — !> leptons. We write it as 

$4 



(2 



"^-^ 1=1,2 



X 5^ {pi - gi - ki)) 5'^ {p2 - 92 - ^2)) 



^ y" dsi ds2 $2 (s, Si, S2) $2 (si, 0, 0) $2 (S2, 0, 0) , 



(177) 



(178) 



(179) 



The form factors Fd^p are computed with s = sh and si^2 — sz to guarantee gauge parameter independence. 
Based on longitudinal polarization dominance (see Appendix |B]) for s — > cxd we introduce longitudinal 
polarization vectors in the numerator of the Z -propagators and obtain 



Y.\iKu+K^)J^J'' =Ek-eL(pi) ^•eL(p2) 



g"^ Re A, 



(180) 



spins 



spins 



A = 



Sz 



1 S - Si - S2 . , ^ 1 SH 

A(S,S1,S2) - o 



COS^0 COS^^ L4 S1S2 -■ "i, -^y - 2 S1S2 

Therefore, taking the residue of all complex poles, in the limit | sh — oo we obtain 



A (s, si, S2) + 4siS2 Fdj 



dsids2 



Sz 



1 



Z'=-i.e+e- i Z<=->n+|j,- 

Ah(s) |2 I Az(si) |2 I AZ(S2) |2 



non-resonant terms. 



(181) 



(182) 



For an Heavy Higgs boson is straightforward to extract the pseudo-observable defining H — ?> Z^Z^ from the 
decay H into two pairs of leptons of different flavor. Interference effects, as well as transverse polarizations, 
will be analyzed in a subsequent paper. 
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